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Abstract 



^ ■ We develop a kinetic theory for point vortices in two-dimensional hydrodynamics. Using 

G ■ standard projection operator technics, we derive a Fokker-Planck equation describing the 

'^ , relaxation of a "test" vortex in a bath of "field" vortices at statistical equilibrium. The 

5 \ relaxation is due to the combined effect of a diffusion and a drift. The drift is shown to be 

O \ responsible for the organization of point vortices at negative temperatures. A description 

that goes beyond the thermal bath approximation is attempted. A new kinetic equation 

is obtained which respects all conservation laws of the point vortex system and satisfies 

Q^ . a H-theorem. Close to equilibrium this equation reduces to the ordinary Fokker-Planck 

equation. 
(N 
l> 

I. INTRODUCTION 

'j^ . It is often useful in two-dimensional turbulence to approximate a continuous vorticity 

field by a cloud of point vortices. The main interest is that such a system is Hamiltonian 
, , ||I|] and can be studied by rather ordinary statistical mechanics. This was first considered 

pj I by Onsager who showed qualitatively the existence of negative temperature states at 

O ■ which vortices cluster. He could therefore explain the occurence of large scales vortices 

(or "supervortices" ) often observed in nature. This was a remarkable anticipation since 
observations were very scarce at that time. His work was pursued by Joyce & Montgomery 
r> \ P] and Lundgren & Pointin who introduced a mean-field approximation and obtained 

c^ ' explicit results for the equilibrium state. They derived in particular a Maxwell-Boltzmann 

statistics for the distribution of point vortices at equilibrium. 

Less is known concerning the relaxation towards equilibium. In fact, the evolution of 
the A^-particle distribution function is governed by a Liouville equation but this equation 
contains too much information to be of practical use. One is more interested by the evolu- 
tion of the one-particle distribution function P(r, t) which gives the probability that a point 
vortex be found in r at time t. In Ref. 0, we have described the relaxation of P{r,t) to- 
wards the Boltzmann distribution in terms of a phenomenological Fokker-Planck equation. 
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In this approach, the vortices have a diffusive motion due to random ffuctuations and they 
experience in addition a systematic drift (Chavanis, 1998) directed along the background 
density gradient. Physically, the drift is the result of a polarization process and its math- 
ematical expression can be determined with a linear response theory ||^. It is found that 
the drift is "attractive" at negative temperatures so the vortices cluster into macrovortices 
in agreement with Onsager's thermodynamical approach. At equilibrium, the drift balances 
the scattering and maintains a non trivial vortex distribution (the Boltzmann distribution) 
providing a dynamical explanation for the persistence of clustering. 

In this paper, we justify our phenomenological model by deriving the Fokker-Planck 
equation directly from the Liouville equation, using projection operator technics [0. These 
methods are standard in statistical mechanics but they are applied here for the first time 
to a system of vortices. We first consider the relaxation of a "test" vortex in a thermal 
bath in which the "field" vortices are in statistical equilibrium. In this approximation, the 
Fokker-Planck equation appears in its usual form with a diffusion term and a drift term. 
The drift coefficient is connected to the diffusion coefficient and to the temperature of the 
bath 1//5 by an Einstein formula. The diffusion coefficient is expressed as a Kubo formula, 
i.e. as the integral of the velocity autocorrelation function. Using an approximation in which 
the vortices are advected by the equilibrium flow, we find that the autocorrelation function 
decays like t~^ for large times 0. This is a slow decay but it ensures the convergence of the 
diffusion coefficient. We also derive non Markovian equations that keep track of memory 
effects. Then, we relax the thermal bath approximation and derive a generalized kinetic 
equation for our vortex system. This integrodifferential equation satisfies all conservation 
laws of the point vortex system and increases the Boltzmann entropy (H-theorem). The 
relaxation towards equilibrium is due to a condition of resonance. If this condition is not 
satisfied, the system can remain frozen in a sort of "metastable" equilibrium. By contrast, 
if the system is sufficiently resonant, it will converge towards the maximum entropy state 
described by the Boltzmann distribution. Close to equilibrium, our generalized kinetic 
equation reduces to the ordinary Fokker-Planck equation. 

The methods developed in this paper are inspired by those introduced in plasma physics 
and stellar dynamics [^. In particular, the systematic drift ^ of a point vortex is the 
counterpart of the dynamical friction [§ experienced by a star in a stellar system. Further 
analogies between 2D vortices and stellar systems are discussed in the paper and in p|-p!4 . 



Other kinetic theories of point vortices have been developed in Refs. [|T5,|T^] in a different 



context. A good review on point vortex dynamics is given by Newton |PT|. 



II. STATISTICAL MECHANICS OF POINT VORTICES 

A. The point vortex model 

In a two-dimensional incompressible fluid, the velocity field u is divergenceless and can 
be written in terms of a stream function %p in the form 

u=-zAV^, (1) 

where z is a unit vector normal to the flow. The stream function is related to the vorticity 
cjz = V A u by the Poisson equation 



00 = -Alp (2) 

obtained by taking the curl of equation (|I]). The impermeabihty condition imphes that ip is 
constant on the boundary and we shall take ^ = by convention. 

We shall consider the situation in which the velocity is created by a collection of A^ point 
vortices of equal circulation 7. In that case, the vorticity field can be expressed as a sum of 
(5-functions in the form 

N 

u{r,t) = Y.^5{v-v,{t)), (3) 

where rj(t) denotes the position of point vortex i at time t. Its velocity is given by 

V, = ^ = -zAV^(r = r,,t), (4) 

where ■?/' is a solution of the Poisson equation (^ with the vorticity field (J^). In an unbounded 
fluid, one has 

1 ^ 
^(r) = -— ^7ln|r-ri|. (5) 

Therefore, the velocity of a point vortex is equal to the sum of the velocities V(j —>■ i) 
produced by the A^ — 1 other vortices, i.e. 

v. = Ev(j-^) (6) 



with 



v(j-.) = --^zA^:^— ^. (7) 

ZTT r,- — r,- r 

The above dynamics can be cast in a Hamiltonian form [jl[: 

dxj _ dH^ dyi_ _ _dl£ . . 

dt dyi ' dt dxi ' 

^=-^E7^1n|r.-r,|, (9) 

where the coordinates (x, y) of the point vortices are canonically conjugate. These equations 
of motion still apply when the fluid is restrained by boundaries, in which case the Hamil- 
tonian (^ is modified so as to allow for vortex images, and may be constructed in terms of 
Green's functions depending on the geometry of the domain. Since H is not explicitly time 
dependant, it is a constant of the motion and it represents the "potential" energy of the 
point vortices (we shall see later on that it also represents the kinetic energy of the flow). 
Therefore, point vortices behave like particles in interaction like electric charges or stars. 
Note, however, that the Hamiltonian (^ does not involve a "kinetic" energy of the vortices 
in the usual sense. This is related to the particular circumstance that a point vortex pro- 
duces a velocity not an acceleration. As a result, an isolated vortex remains at rest contrary 
to a material particle which has a rectilinear motion due to its inertia. Point vortices form 
therefore a very peculiar Hamiltonian system. 



B. The microcanonical approach of Onsager 

The statistical mechanics of point vortices was first considered by Onsager ^ who showed 
the existence of negative temperatures at which point vortices cluster into "supervortices" . 
Let us briefiy recall his argumentation. 

Consider a liquid enclosed by a boundary, so that the vortices are confined to an area A. 
Since the coordinates (x, y) of the point vortices are canonically conjugate, the phase space 
coincides with the configuration space and is finite: 

/ dxidyi...dx]\fdyN = i / dxdy) = A^ . (10) 

This striking property contrasts with most classical Hamiltonian systems considered in sta- 
tistical mechanics which have unbounded phase spaces due to the presence of a kinetic term 
in the Hamiltonian. 

As is usual in the microcanonical description of a system of A^ particles, we introduce 
the density of state 

g{E) = J dxidyi...dxNdyNS(E - H{xi,yi,...,XN,yN)] (H) 

which gives the phase space volume per unit interaction energy E. The phase space volume 
which corresponds to energies H less than a given value E can be written 

$(E) = / g{E)dE. (12) 

It increases monotonically from zero to A^ when E goes from Emin to +00. Therefore, 
g{E) = d^{E)/dE will have a maximum value at some E = Em, say, before decreasing to 
zero when E ^ +00. 

In the microcanonical ensemble, the entropy and the temperature are defined by 

S = \ng{E) P = ^ = ^- (13) 

For E > Em, S{E) is a decreasing function of energy and consequently the temperature 
is negative. Now, high energy states E 3> Em are clearly those in which the vortices 
are crowded as close together as possible. For energies only slighlty greater than Em, the 
concentration will not be so dramatic but there will be a tendency for the vortices to group 
themselves together on a macroscopic scale and form "clusters" or "supervortices". By 
contrast, for E < Em, the temperature is positive and the vortices have the tendency to 
accumulate on the boundary of the domain in order to decrease their energy. For a system 
with positive and negative vortices, the negative temperature states, achieved for relatively 
high energies, consist of two large counter-rotating vortices physically well separated in the 
box. On the contrary when E -^ —00, the temperature is positive and vortices of opposite 
circulation tend to pair off. 



C. The mean field approximation 

It is easy to show that the exact distribution of point vortices expressed in terms of 
5-functions 

N 
1=1 

is solution of the Euler equation 

" + UexVUex = 0, (15) 



dt 

where u^x is the exact velocity field determined by equations ([^) (0) and (0). This is proved 
as follows. Taking the derivative of ([l^ with respect to time, we obtain 



^ = -E7V5(r-r,(t))V,. (16) 

Using Vj = Uex{ri{t),t), we can rewrite the foregoing equation in the form 

^ = -VE75(r-r.(t))u..(r,t). (17) 

Since the velocity is divergenceless, we obtain 

^(r,t) = -Ue,.(r,t)VE75(r-r.(i)) = -u,xVuj,x{r,t). (18) 

Therefore, the Euler equation (|15D with (|l]) (^ and (|^) contains exactly the same information 
as the Hamiltonian system (|§)(^). 

This description in terms of 5- functions, while being technically correct, is useless for 
practical purposes, because it requires the knowledge of the exact trajectories of the point 
vortices for an arbitrary initial condition (or the solution of the Euler equation (|15|)). When 
N is large, this task is impossibly difficult. Therefore, instead of the exact vorticity field 
expressed in terms of delta functions, one is more interested by functions which are smooth. 
For that reason, we introduce a smooth vorticity field {uj){r,t) which is proportional to 
the average number of vortices contained in the cell (r, r + dr) at time t. This mean field 
description, which ignores the granularities of the sytem, requires that it is possible to divide 
the domain in a large number of cells in such a way that each cell is (a) large enough to 
contain a macroscopic number of point vortices but (b) small enough for all the particles in 
the cell can be assumed to possess the same average characteristic of the cell. 

In this mean field approximation, the Hamiltonian (^ is changed into 

E = -— f {iu){r){uj){r')\n\v - r'\dhdh'. (19) 

In writing this expression we have not taken into account the constraint j ^ i appearing in 
(P). Really, in (|19D the integration extends over the point r = r' so that ([T9|) contains self- 
energy terms which become infinitely large for point vortices. As will soon become apparent, 



this meanfield approximation implies that the energy E is positive, a property which is not 
necessarily shared by the Hamiltonian @. 

Using equation (|^), adequatly generalized to account for a continuous distribution of 
vortices, our expression (|19D for E can be rewritten 

E = ^J{uj)i,dh. (20) 



Introducing explicitly the Poisson equation in ( poD and integrating by parts, one has succes- 
sively 

E = lJ ^{-A^)dh = \ JiV^fdh = J ^-^d\ (21) 

where (u) is the smooth velocity field (the second equality is obtained by a part integration 
with the condition ?/; = on the boundary). Therefore, E can be interpreted either as the 
potential energy of interaction between vortices (see equation (pOl)) or as the kinetic energy 
of the flow (see equation (pl|)). 

D. The mean field equilibrium 

We now wish to determine the distribution of vortices at equilibrium following a statis- 
tical mechanics approach [Q. To that purpose, using Boltzmann procedure, we divide the 
macrocells (r, r + dr) into a large number of microcells and enumerate the number of "mi- 
crostates" which correspond to the same "macroscopic" configuration of the system. The 
logarithm of this number defines the entropy. In the mean field approximation, this leads 
to the classical formula 

S = ~N f P(r) \nP{r)dh, (22) 

where -P(r) is the density probability that a point vortex be in the surface element centered 
on r. The average vorticity in r is related to this probability density by 

{uj){r) = N-fP{r). (23) 

At equilibrium, the system is in the most probable macroscopic state, i.e. the state that is 
the most represented at the microscopic level. This optimal state is obtained by maximizing 



the Boltzmann entropy (^) at fixed energy (poD and vortex number N, or total circulation 



T = N-f= {uj)dh. (24) 

Writing the variational principle in the form: 

6S - I36E - a6T = 0, (25) 

where j3 and a are Lagrange multipliers, we find that the maximum entropy state corresponds 
to the Boltzmann distribution 1311 : 



(uj) = Ae-f^^^ (26) 

with inverse temperature (3. We can account for the conservation of angular momentum 
L = J{uj)r'^(Pr (in a circular domain) and impulse P = ^ {uj)yd?Y (in a channel) by intro- 
ducing appropriate Lagrange multipliers Vt and U for each of these constraints. In that case 
equation (^) remains valid provided that we replace the streamfunction ip by the relative 
streamf unction if)' = il) + ^r'^ — Uy. This more general situation has been considered in, e.g., 



Ref. [|T8l to describe rotating or translating dipoles. 

Substituting the Boltzmann relation between (ui) and ip in the Poisson equation (|^), we 
obtain a differential equation for the streamfunction: 

- Alp = Ae"^'^'^ {point vortices) (27) 

which determines the equilibrium distribution of vortices. In the case of stellar systems and 



electric charges, the corresponding Boltzmann-Poisson equation has the form [|T9 



A$ = 47rGv4e-^'™* {stellar systems) (28) 

A 
- A<l> = — e"^'"^* {electric charges) (29) 

Co 

where $ denotes successively the gravitational and the electrostatic potential. For these 
systems, (3' > since the temperature is a measure of the kinetic energy. By contrast, for 
point vortices there is no kinetic term in the Hamiltonian (^ and the temperature can be 
either positive or negative. When f3 < 0, equation (^) is similar in structure to equation 
(p8|). The vortices tend to attract each other, like stars in a galaxy, and form "clusters" or 
"supervortices" . The density profile determined by (|27| ) or (|28| ) is a decreasing function of the 
distance. When P > 0, equation (pTl) is similar in structure to equation (|29|) . The vortices 
tend to repell each other, like electric charges, and accumulate at the boundary. The density 
profile determined by (^) or (p9|) is an increasing function of the distance. Therefore, the 
formal analogy between 2D vortices and stellar systems is intimately related to the existence 
of negative temperatures in 2D turbulence. However, the physical mechanism by which 
vortices and stellar systems achieve equilibrium is different. Whereas the organization of 
stars is relatively clear because of the attractive nature of gravity, the organization of point 
vortices at negative temperatures is much less intuitive. In the following section, we shall 
give a physical interpretation of this phenomenon in terms of a "systematic drift" . 



III. ELEMENTARY DERIVATION OF THE SYSTEMATIC DRIFT 

A. Analogy with Brownian motion: the necessity of the drift 

We shall first show the necessity of this drift by using an analogy with Brownian theory. 
The starting point of this analogy is to realize that the velocity of a point vortex can be 
decomposed in two terms: a smoothly varying function of position and time (V)(r, t) and 
a function V{t) taking into account the "granularity" of the system and undergoing strong 
discontinuities. The total velocity of a point vortex can therefore be written: 



V=(V)(r,t)+V(t). (30) 

The velocity (V)(r, t) reflects the influence of the system as a whole and is generated by the 
mean vorticity (u) (r, t) according to the Biot & Savart formula: 

(V)(r,t) = -^zA /^;-^M(r',t)dV. (31) 

Zir J |r — r|^ 

The fluctuation V(t) arises from the difference between the exact distribution of the point 
vortices ujexact{^,t) and their "smoothed-out" distribution {uj){r,t). It is on account of these 
fluctuations that the velocity of the test vortex will depart from its mean field value (V). 
The velocity fluctuation V, of order ^ (where d ~ n~^/^ is the inter- vortex distance), is much 
smaller than the average velocity (V), of order wfR (where R is the domain size), but this 
term has a cumulative effect which gives rise to a process of diffusion. It makes sense therefore 
to introduce a stochastic description of the vortex motion like for colloidal suspensions in a 
liquid [^ or stars in globular clusters [§]. However, contrary to the ideal Brownian motion, 
point vortex systems have relatively long correlation times. This makes the study much 
more complicated than usual and the technical study of section |rV| is required. In order to 
gain some physical insight in the problem, we shall ignore this difficulty for the moment and 
describe the system by traditional stochastic processes. 

According to equation (0), we would naively expect that the evolution of the density 
probability P(r,t) be governed by a diffusion equation of the form 

dP 

— + (V)VP = DAR (32) 

This would in fact be the case for a passive particle having no retroaction on the vortices 
or when the distribution of vortices is uniform like in [^ |jl3[| |2^. However, this diffusion 
equation cannot be valid when the system is inhomogeneous. A first apparent reason is that 
equation ( |32D does not converge towards the Boltzmann distribution ( p6D when t — ^ +oo. 
Another related difficulty is that equation (^) does not conserve energy. It seems therefore 
that a term is missing to act against the diffusion. 

These problems are similar to those encountered in Brownian theory or for stellar sys- 
tems. They have traditionally be solved by introducing a dynamical friction in order to 
compensate for the effect of diffusion. The occurence of this frictional force is a manifesta- 
tion of the "fluctuation-dissipation" theorem in statistical mechanics. In the present context, 
the dynamical friction is replaced by a systematic drift of the vortices. We must therefore 
rewrite the decomposition (|30|) in the form 

V=(V)-eV^+V(t), (33) 

where ^ is the drift coefficient. In section [1IIB| , we shall give a physical justification for the 
existence of the drift in terms of a polarization process and in section |^ we shall derive 
this term directly from the Liouville equation by using projection operator technics. The 
importance of this drift was first pointed out by Chavanis (1998) using a thermal bath 
approximation and a linear response theory. The drift term must be calculated by resorting 
to relatively elaborate technics but it is remarkable that a general relationship between ^ 
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and D can be obtained without beeing required to analyze at any point the details of the 
"sub dynamics" . 

According to equation (p3D, the equation of motion for a point vortex can be written in 
the form 

^ = (V)-eVV^+V(t). (34) 

Since the velocity V(t) undergoes strong discontinuities, the trajectory r(t) of the point 
vortex is not difFerentiable. Therefore, equation (|3^) must be viewed as a stochastic equation 
analogous to the Langevin equation in the ordinary Brownian theory. Let At be an interval 
of time long compared to the fluctuation time but short at the scale on which the physical 
parameters change appreciably. The variation in the position Ar of the particle during At 
is given by 

Ar = ( V) At - ^ V^ At + B ( At) , (35) 

where 

B(At) = / V{t')dt'. (36) 

Each fluctuation V produces a small displacement 5v but the repeated action of these fluc- 
tuations produces a net displacement of the same order as the drift ^Vip. To determine the 
probability iL'[B(At)] that the fluctuations produce a displacement B(At) during the time 
interval At, we first divide the interval (t, t + At) into a succession of discrete increments 
in position and observe that B(At) is a sum of N random variables T(Vj)Vi where T(V) 
characterizes the typical duration of the velocity fluctuation V. This is a problem of random 
walks where B(At) represents the distance reached after N steps. For large A^'s, the Central 
Limit Theorem leads to a Gaussian transition probability: 

with a diffusion coefficient 

D = \{T{V)V'). (38) 

We now assume that the motion of a point vortex can be idealized by a Markov process, 
i.e. the probability at time t + At depends on the probability at time t but not at earlier 
times. As indicated previously, this approximation is not completely correct in the case of 
point vortices which have long correlation times. However, relaxing this hypothesis would 
involve more intricate equations (see section [TVD and we shall ignore this difficulty for the 
moment. We write therefore: 

P(r, t + At) = j P{y - Ar, t)w(r - Ar| Ar)d2( Ar), (39) 

where w{v — Ar|Ar) is the probability for a point vortex located in r — Ar to suffer an 
increment of position Ar during At. Expanding P(r,t + At), P(r — Ar, t) and i(7(r — Ar|Ar) 
in the form of Taylor series, we arrive at the Fokker-Planck equation in its general form: 



where 



fA*=-E|-(^(A''.»4E5|^(^(A''.A'',». (40) 



{An) = f Ariw{r\Ar)d^{Ar) (41) 



{AriAvj) = f AriArjw{v\Ar)d'^{Ar). (42) 

According to equation (P7|), the transition probabihty from r to r + Ar is given by 

w(r Ar) = —expl ^^ ' , ^^ — ^ L (43) 

^ ' ' inDAt ^l iDAt J ^ ^ 

With (^3]), the moments (^) (|4^ ) can be easily evaluated yielding 

(Ar) = ((V) - eV^)At, (Ar^Ar^) = 2DAt6ij. (44) 

Substituting these results in the general Fokker-Planck equation (|40|), we find that: 



9P 

— + (V)VP = V(Z)VP + ePVV^). (45) 

We had previously introduced this equation in Ref. p using phenomenological arguments. 
The physical interpretation of each term is straightforward. The left hand side (which can 
be written dP/dt) is an advection term due to the smooth mean field velocity (V) . The right 
hand side can be written as the divergence of a current —V.J and is the sum of two terms: 
the first term is a diffusion due to the erratic motion of the vortices caused by the fluctuations 
V; the second term accounts for the systematic drift of the vortices. At equilibrium, the drift 
precisely balances random scatterings and the distribution ( p6D is settled. More precisely, 
the condition that the Maxwell-Boltzmann statistics ( PB| ) satisfles equation (^) identically 
requires that D and ^ be related according to the relation 

e = D(3^ (46) 

which is a generalization of the Einstein formula to the case of point vortices. A more rigorous 



justification of this relation will be given in section |^ where the diffusion coefficient and 
the drift term are calculated explicitly. 

B. Systematic drift: the result of a polarization process 

According to the previous discussion, the relaxation of a point vortex towards statistical 
equilibrium can be described by a Fokker-Planck equation 

dP 

— + (V) VP = V{D{VP + /?7PVV')) (47) 
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involving a diffusion term —DVP and a drift term 

{y)dr.ft = -D/3-fV^. (48) 

The drift is normal to the mean field velocity (V) = —zAVip of the vortices and its direction, 
depending on the sign of /?, has important physical implications. To fix the ideas, let us 
assume that all point vortices have positive circulation (the opposite case leads to the same 
conclusion). Due to the meanfield velocity, a particular point vortex rotates anticlockwise. 
At negative temperatures, the drift is directed to its left and the vortex is attracted to the 
center of the domain. On the contrary, at positive temperatures, the drift is directed to its 
right and the vortex is rejected against the boundary. Therefore, the effect of the drift is 
consistent with Onsager thermodynamical approach and it provides, in addition, a physical 
mechanism to understand the organization of point vortices at negative temperatures. For 
/? = 0, the medium is homogeneous and there is no drift. Equation (|T7|) reduces to a pure 



diffusion equation like in Ref. pi| , p!3| , |22| . Therefore, the drift occurs only in the presence of 
a background shear. 

In stellar systems, the relaxation of the distribution function /(r, v, t) is usually described 
by the Kramers- Chandrasekhar equation: 



which is a particular Fokker-Planck equation with a structure analogous to equation 

In this analogy, the dynamical friction experienced by a star as a result of close encounters 

(Chandrasekhar IQ): 

(F) friction = -Dpmv (50) 



is the counterpart of the systematic drift (^q ) experienced by a point vortex in two dimen- 
sional turbulence. The dynamical friction can be viewed as the drag exerted on a test star 
by the wake it induces in the field stars, like in a polarization process. We can use a simi- 
lar approach to understand the origin of the drift. Let us consider a collection of A^ point 
vortices at statistical equilibrium with inverse temperature /3. When (3 < 0, the density of 
these "field" vortices decreases from the center to the periphery of the domain. A "test" 
vortex moving through this medium locally modifies the vorticity distribution and produces 
a polarization cloud which amounts typically to a rotation of the surrounding vortices. This 
creates an excess of density behind it an a deficit of density in front of it. Therefore, the 
retroaction of the field vortices leads to a drift of the test vortex directed inward. We reach 
the opposite conclusion if the temperature is positive. When the system is homogeneous 
(/3 = in a domain with no special symmetry), the polarization cloud induced by the test 
vortex has no effect and the drift cancels out. Therefore, a homogeneous system of point 



vortices remains homogeneous [21,13,22 



C. The Maximum Entropy Production Principle 

In section ^V\, we shall derive the Fokker-Planck equation ( P7| ) directly from the Liouville 
equation. However, we want to show first that the general structure of this equation can be 
understood from relatively simple thermodynamical arguments. 

11 



Let us rewrite equation (^) in the form 

^ + (u)V(a;) = -V.J, (51) 

where J is an unknown diffusion current. This equation conserves the circulation (p4D 
provided that J.fi = on the domain boundary (fi is a unit vector normal to the boundary). 
The problem at hands consists in determining the expression for J. Its exact expression 
depends on microscopic processes and is therefore difficult to capture. However, it is easy 
to write down some macroscopic constraints that it must satisfy. These constraints are 
provided by the first and second principles of thermodynamics, namely the conservation of 
energy and the increase of entropy. We shall find that these constraints are very stringent 
and determine completely the structure of the diffusion current. 



Taking the time derivative of equations (^) (p2D and substituting for (^) we obtain the 
constraints 



E= J.Vipd'r = (52) 

and 

S= / J. V \n{uj)dh > 0. (53) 

7 J 

We shall now introduce an optimization procedure known as the Maximum Entropy Pro- 
duction Principle (M.E.P.P.). This principle was introduced initially in the context of 2D 
turbulence by Robert & Sommeria ^3\ but its domain of applicability is very general and 



concerns, for example, the case of stellar systems P,|2^. This principle states that "out of 



equilibrium, the system evolves so as to maximize its rate of entropy production S while 
accounting for all the constraints imposed by the dynamics, in particular the conservation 
of energy E = 0". There is no precise justification for this principle and it is important 
therefore to confront the M.E.P.P. with more rigorous methods, like those of section |V|, to 
determine its domain of validity. In any case, the M.E.P.P. can be considered as a conve- 
nient tool to build relaxation equations which are mathematically well-behaved and which 
can serve as numerical algorithms to calculate maximum entropy states. 

We seek therefore the optimal diffusion current J which maximizes the rate of entropy 
production S at fixed energy. In order to avoid the unphysical solution |J| — > +oo with 
J.Vip = 0, we impose the additional constraint 

j2<C(r,t), (54) 

where C is an upper bound which must exist but is not known. The solution of the opti- 
mization problem is 

J = -D{V{u) + /3(t)7(cu)V^), (55) 

where /3(t) and D{r,t) are Lagrange multipliers associated with the constraints (|52D and 
(^). When substituted in equation (0), we obtain: 

^ + (V) V(a;) = V(D(V(c.) + /?(t)7(a;) W^)). (56) 
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This equation has the same form as the Fokker-Planck equation (|47|). Here, the diffusion 
term arises from the variations of entropy 6S and the drift term is necessary to conserve 
energy. Note that the Einstein formula ( ^B]) is automatically satisfied by this variational 



approach. 

The time evolution of the inverse temperature /3(t) is determined by the conservation of 
energy. Substituting the diffusion current (E^ in the constraint (|52D we find 



We can also check that the entropy monotonically increases during the relaxation provided 
that D > 0. Indeed, using (|55D(|55D and ( [5^ ) we can easily establish that 



'-Id^/''^°- («^' 



At equilibrium, J = and we recover the Boltzmann distribution (pGI). 

Note that the optimal current (^) can be written J = x^a where a = \n{uj) + P'yip is 
a "generalized potential" which is uniform at equilibrium. Therefore, the M.E.P.P. can be 
viewed as a variational formulation of the linear thermodynamics of Onsager which relates 
the diffusion currents to the gradients of generalized potentials. However, the M.E.P.P. gives 
a more elegant approach to the problem and, conceding a real importance to the constraints, 
it is easier to implement in more comphcated situations |2^J^,^ . In addition, it shows that 
the structure of the relaxation is determined by purely thermodynamical arguments. All 
explicit reference to the subdynamics is encapsulated in the diffusion coefficient which is 
left unspecified (it appears as a Lagrange multiplier related to an unknown bound on the 
diffusion current). It must be therefore calculated with a more microscopic model like the 
one of section |^. 

IV. RELAXATION OF A POINT VORTEX IN A THERMAL BATH 

A. The Liouville equation 

Let us consider a collection of A^ + 1 point vortices with identical circulation 7. We 
select one of these vortices, for example point vortex 0, and call it the "test vortex". The 
other vortices 1, ...,N will be refered to as the "field vortices". Let yu(r,ri, ...,riy,t) denote 
the A^ + 1 particle distribution of the system, i.e /i(r,ri, ...,riy,t)cPr(fri...cPrN represents 
the probability that point vortex be in the cell (r, r + dr), point vortex 1 in the cell 
(ri, ri + (iri)... and point vortex A^ in the cell (r^v, rjv + rfrjv) at time t. The (A^ + l)-particle 
distribution function fi(t) satisfies the Liouville equation 

| + f;v.|^=0, (59) 

dt ^ dvi 

where Vj is the velocity of vortex i produced by the other vortices according to equations 
(i) (0)- We also introduce the one- and A^-particle distribution functions defined by 
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N 

Pir,t)= J fii{rk},t)lldhk, (60) 



fisys{ri, ..., Tn, t) = J f^ii^k}, t)d^r. (61) 

The physical picture that we have in mind is that the test vortex evolves in a "bath" of field 
vortices. Therefore, we rewrite the distribution function fi in the suggestive form 

/x(r, ri, ..., Tat, t) = P(r, t)fisys{ri, ..., tn, t) + ;U/(r, ri, ..., yn, t), (62) 

where /i/ refiects the effect of correlations between the test vortex and the field vortices. 
Physically, this term accounts for the polarization process described qualitatively in section 

EH 

The Liouville equation (|59D provides the correct starting point for the analysis of the 
dynamics of our vortex system. However, when N is large, this equation contains much 
more information than one can interpret. Consequently, what one would like to do is to 
describe the system in some average sense by a one-particle distribution function. In the 
previous sections we have indeed derived heuristically some differential equation satisfied by 
this distribution function on the basis of stochastic arguments. We shall now discuss the 
connection of such heuristic theories with a more microscopic description of the system. 



B. The projection operator formalism 

Our first objective is to derive some exact kinetic equations satisfied by P(r, t) and 
fisysi'i'i, ..., iLN, t)- This can be achieved by using the projection operator formalism developed 
by Willis & Picard . This formalism was also used by Kandrup |^ in the context of stellar 
dynamics to derive a generalized Landau equation describing the time evolution of the 
distribution function of stars in an inhomogeneous medium. We shall just recall the main 
steps of the theory. More details can be found in the original paper of Willis & Picard 
and in Kandrup [^. To have similar notations, we set x = {r} and y = {ri, ...,rN}. Then, 
equation (p2| ) can be put in the form 

Kx, y, t) = fXR{x, y, t) + fj,i{x, y, t) (63) 

with 

Ij,r{x, y, t) = /(x, t)g{y, t), (64) 

where we have written f{x,t) = -P(r, t) and g{y,t) = fisys{ri, ...,rN,t). The Liouville equa- 
tion is also cast in the form 

du 

— = -iLfi = -i{Lo + Lsys + L')^, (65) 

where Lq and Lgys act respectively only on the variables x and y, whereas the interaction 
Liouvillian L' acts upon both x and y (the complex number i is here purely formal and has 
been introduced only to have the same notations as Ref. PJ^l)- 

14 



Following Willis & Picard, we introduce the time-dependant projection operator: 

P{x, y, t) = g{y, t) J dy + f{x, t) J dx - f{x, t)g{y, t) J d^ J dy. (66) 

We can easily check that: 

P(x, y, t)fi{x, y, t) = fiR{x, y, t), (67) 

[1 - P(x, y, t)]/i(x, y, t) = ni{x, y, t). (68) 

We also verify that P is a projection in the sense that P^(t) = P(t). Applying P and 1 — P 
on the Liouville equation (^), we obtain the coupled equations 

dtHnix, y, t) = -iPLfiR - iPLfxj (69) 

and 

dtfiiix, y, t) = -i{l - P)Lfin - ^(1 - P)Lfii. (70) 

These equations should be compared with equations (8) (9) that appear in the quasilinear 
theory of 2D turbulence |2^. In the present context, equations ( |69D ( ffOD describe the separa- 



tion between a "macrodynamics" and a "subdynamics" . In the quasilinear theory, equations 
(8) (9) describe the evolution of the "coarse-grained" and "fine-grained" components of the 
vorticity. 

Introducing the Greenian 

git,t') = expl-t £dt"[l - P{t")]LV (71) 

we can immediately write down a formal solution of equation (^), namely: 

mix, y,t) = - f dt'Git, t')i[l - Pit')]Lf^Rix, y, t'), (72) 

Jo 

where we have assumed that initially the particles are uncorrelated so that ^j{x,y, 0) = 0. 
Substituting for fij{x,y,t) from equation ([7^ ) in equation (p^), we obtain 

dtfXR{x,y,t) = -iPL^iR- f dt'P{t)LG{t,t')[l- P{t')]L^iR{x,y,t'). (73) 

The integration over y will yield an equation describing the evolution of /. Using some 
mathematical properties of the projection operator (|66D, the final result can be put in the 
nice symmetrical form given by Willis & Picard [||: 

dtf{x, t) + iLof + i{L')sysf = - J^ dt' J dyAtL'Git, t')AeL'g{y, t')f{x, t'), (74) 

where the notations stand for 
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{L'),y, = ldy'L'{x,y')g{y',t), (75) 

{L')o = ldx'L'ix',y)fix',t), (76) 

AtL' = L'- {L')sys - {L')o. (77) 

Similarly, after integrating over x we find the equation satisfied hj g 0: 

dtg{y, t) + iLsysg - t{L')ig = - j dt' j dxAtL'g{t, t')At>L'g{y, t')f{x, t'). (78) 

C. Application to the point vortex system 

The previous theory is completely general and we now consider its application to a system 
of point vortices. Let us first rewrite the Liouville equation ([59| ) in a form that separates 
the contribution of the test vortex from the contribution of the field vortices: 

|,Jv,-.o)|.Ev(o..,| 

+ E E V(j-z)|^ = 0. (79) 



The different operators that arise in the decomposition (65) are 



iLo = 0, (80) 

N Q 

iLsys = Y.ll^U^^)Trr. (81) 

d _.„ .. d 



■i' = g{v(,^0)- + V(0^.)^}. (82) 

The mean-field velocity created by the field vortex i on the test vortex is denoted by 

(V(z ^ 0)) = I P(r„ t) V(2 ^ 0)d\. (83) 

Similarly, 

(V(0 ^i))= j P(r, t)V(0 -^ i)dh (84) 

denotes the mean field velocity created by the test vortex on the field vortex i. Finally, the 
total mean field velocity experienced by the test vortex is given by 
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N N 

(V) = E(V(^ - 0)) = 5:(V(1 -. 0)) = iV(V(l -. 0)), (85) 

where the second equahty follows from the identity of the point vortices. 

We are now ready to evaluate the quantities (^) (|T^) ([77|). After straightforward inte- 
gration by parts, we find successively 

f\ N a 

^{L')sy. = {y)g^ = T.{Vi^^O))-, (86) 

-/V o 

^(^% = E(V(0-z))^, (87) 

N f "^ Pi ^ I" "^ fi 

zA,L' = E V(z ^ 0) - (V(. ^ 0)) - + E V(0 - .) - (V(0 -. ^)) — . (88) 

Introducing the velocity fluctuations: 

V{i ^ 0) = V(z ^ 0) - (V(z ^ 0)), (89) 

V(0^z)=V(0^i)-(V(0^i)), (90) 

we can rewrite our expression for AfL' in the form 

N Cj N a 

^A,L' = EV(z - 0)- + EV(0 - z) — . (91) 

Substituting these results in equation (0), we obtain the following kinetic equation for 
the one-particule distribution function of a vortex system: 



r)P riP rt r ^ ^ r) 



x^(t,t') EV(j-0)|^ + EV(0^j)t^ P(r,t')/^.,,({r,},t') (92) 

or, alternatively, 

dP „..dP a r' - " " " 



^^ ^' •^'^ -' fc=i i=ij=i 



where the Greek indices refer to the components of V in a fixed system of coordinates. We 
can note that equation (P3|) already shares some analogies with the Fokker-Planck equation 



of section III. Indeed, the first term on the r.h.s. corresponds to a diffusion and the second 



term to a drift. For a passive particle V {0 —^ j) = and the drift cancels out, as expected. 
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D. The thermal bath approximation 

Equation (^) is an exact differential equation for P{r,t). However, tliis equation is not 
directly soluble since tlie unkown function fisys{{''^k},t) is given by an equation of tlie form 
([78|) depending in turn on P(r, t). We tlierefore liave to solve the coupled system (|7^ -(|78D. 
This system bears exactly the same information as the initial Liouville equation (^) and, 
without further simplification, is untractable. 

To reduce the complexity of the problem, we shall implement a "thermal bath approxi- 
mation" . We assume that the field vortices are in statistical equilibrium with inverse temper- 
ature Peg- Therefore, the iV-particle distribution function ^sys{{T'k}) can be approximated 
by a product of A^ one-particle distribution functions Peg given by the Maxwell-Boltzmann 
statistics (^). In other words, we make the approximation 

llsysin, ■■■, Tat, t) ~ fJ.eq{ri, ■■■, ^n) (94) 

with 

Af N 

/ie,(ri,...,r^) = n^e,(rfe) = J] ^-^e'^-^'^- ('•'=), (95) 

fc=i fc=i 

where ipeq is solution of the Poisson equation (g) with the equihbrium vorticity {uj)eq = 
N'jPeq. Substituting explicitly for the Boltzmann distribution from equation (|95D in equation 
(PB|), we obtain 

dP dP d rt r ^ ^ ^ 

xe(«,« - r) (yU - 0) A - A,7V''(0 ^ j)^)P(r, i - t) H ft,(r»), (96) 

where {.)eq denotes the average with respect to the equilibriun distribution Peg. Explicating 
the action of the Greenian, we can rewrite our equation ( p6l) in the form 

N N N 



^ + (V)e.¥ = ^ Hr / n rf^r, ^ ^ ^^(z - 0, t) 






x(v^(j-0,t-r)— -/3.,7V^(0-j,t-r)^(r,(t-r))) 

xP(r(t - r), t - r) n ^e,(rfe(t - r)), (97) 



fe=i 



where the retarded velocity V{j ^ i,t — t) must be viewed as an explicit function of time. 
More precisely V{j ^ i,t — t) is a shorthand notation for V{rj{t — t) ^ Yi{t — r)), where 
Yi{t — t) denotes the position at time t — r of the i-th point vortex located in Tj = rj(t) 
at time t. The trajectories of the point vortices between t — t and t are determined by 
the complicated Greenian Q(t,t — r) defined in ([711). We need therefore to solve the exact 
Kirchhoff-Hamilton equations of motion (Rl). In fact, to a good approximation, we can 



consider that the point vortices are purely advected by the equihbrium mean field velocity 
(V)eg. Indeed, when N -^ oo, we have already indicated that the velocity fiuctuation 
V is much smaller than the mean field velocity {V)eq- Therefore, we can replace the exact 
Greenian ^ by a smoother Greenian {G)eq which would be obtained if the point vortices were 
moving in the velocity field created by the equilibrium distribution function fieq- Formally, 
this Greenian is constructed with the averaged Liouville operator {L)eq = J2iLo0^i)eq-^- 

In this approximation, the correlations involving two different vortex pairs vanish and the 
equation can be simplified considerably. Using the results of Appendix ^ (see in particular 
equation ( |B11| )), we find 



xP(r(t-r),t-r)Peg(ri), (98) 

where we have used Peqijiit — T)) = Peq{ri) since Peq = fi^peq) is constant along a streamline 
and the particles are assumed to follow the streamlines in a first approximation. Since the 
vortices are identical we also have 

dP ,_, dP _ d 



x(K^(1^0,t-r)A_;3e,7nO-l,t-r)^(r,(t-r))) 

xP(r(t-r),t-r)Pe,(ri). (99) 



Noting that the integral is dominated by the divergence of the product V^V^ when ri -^ r, 
we can make a "local approximation" and replace -^^(ri) and Peg(ri) by their values taken 
in r. For the same reason, we can neglect vortex images and replace the Kernel V(0 -^ 1) 
by its singular part (^ satisfying V(0 — i> 1) = — V(0 -^ 1). With these approximations, the 
kinetic equation (PO) takes the form: 



dt 

X 



-r 



8P f)P f) r^- r 

■|^(r(t - r), t - r) + (3^,^P{r{t -r),t- r)^(r(t - r))) . (100) 



E. The Fokker-Planck equation 

1. Unidirectionnal flow 

We will now see how the previous equation can be simplified for particular equilibrium 
fiows. We shall first consider the case of an unidirectional fiow {V)eq = Kg(2/)x produced by 
a vorticity distribution {uj)eq{y)- If we restrict ourselves to solutions of the form P = P{y, t), 
the kinetic equation ( |100| ) becomes 
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dP 



d 



Or u t'- f 

-Q^ = ^Q-ldrJ d^v,P,,{y)Vy{l ^ 0, t)^^! ^ 0, t - r) 



X (^{y, t-T)+ f3,,^P{y, t - T)^{y) 



where we have used y{t — T) = y{t) = y since the point vortices follow the streamlines of the 
equilibrium flow. This equation can be rewritten in the form 

at 

where 



aWo'^^^^)fe^^'^-^^ + ^-^^^^'^-^^^^ 



(102) 



C{r) = C™(r) =nJ dhiVy{l -^ 0, t)Vy{l ^ 0, t - r)P,g{y) 
is the velocity autocorrelation function. In Appendix P 1| , it is found that: 



C{t) 



"StT 



In AT 



l + iS2(y)r2 ^^ 



PeM, 



where 



^iy) = -^{v)eM 



(103) 



(104) 



(105) 



is the local shear of the flow, equal here to the vorticity. Note that the velocity autocorre- 
lation function decays like r~^ for r —>■ +oo P). 

Equation (|102| ) is a non Markovian equation since the probability P{y, t) at time t de- 
pends on the probability P{y, t — t) at earlier times through an integration over r. Accord- 
ingly, the present study which explicitly takes into account memory effects is more general 
than the stochastic model presented in section pi A| . However, if we implement a Markov 
approximation and replace P{y, t — r) by P{y, t), we recover the Fokker-Planck equation of 



section 111: 



dP 

'dt 



d_ 
dy 



V oy oy 



with a diffusion coefficient 

D = Dyy = N 

and a drift term 



j^ drj dh,Vy{l ^ 0, t)l^^(l -^ 0,t - T)P,g{y) 



Vy 



■iVy) 



drift 



l3eqlD 



dy ' 



(106) 



(107) 



(108) 



The drift coefficient is given by an Einstein relation as expected from the general consider- 
ations of section [111 A| . The diffusion coefficient is expressed as a Kubo formula, i.e. as the 
integral of the velocity correlation function (see Appendix ^ . In Appendix |C1| it is found 
that (see also Ref. (O): 



D = -N-f^ 



1 



|S(l/)| 



\nNP,,iy). 



(109) 



The reason for the logarithmic divergence is explained in Appendix ^ and in Ref. [|T3 
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2. Axisymmetrical flow 



We now consider the case of an axisymmetrical equilibrium flow such that {V)eq = 
{V)eq{r)e0. This flow is generated by an equilibrium vorticity field {uj)eq{r). If we restrict 
ourselves to solutions of the form P = P{r,t), the kinetic equation (|100|) simplifies in 



dP 



'L_ = N-^r I dT f dhiPegir)^^'^! ^ 0,t)V'''^'-^\l -^ 0,t 



/dP 
\ dr 



r,t-T)+ f3eqlP{r, t - T 



d^, 



eq 



dr 



(110) 



where we have used r{t — t) = r{t) = r since the point vortices follow the streamlines of the 
equilibrium flow. Furthermore, \/^'(*) and \/'"(*~'^) denote the radial components of V(l — »• 0) 
at times t and t — r in a polar system of coordinates (ej.(t), e6i(t)) moving with the test 
vortex. Equation (TH,) can be rewritten 



where C(r) is the velocity autocorrelation function 

C[t) = N fdhiV^^'\l -^ 0,t)V"(*-")(l ^ 0,t -r)Peq(r). 
In Appendix |C2| it is found that 






l + iS2(r)r2 "^ 



Peq{r), 



where 

dr\ r 

is the local shear of the flow. 

If we ignore memory effects, we obtain the Fokker-Planck equation: 

dP _1 d 

dt r dr 



rD\-:—+ PeqlP 



V dr 



dr 



with a diffusion coefficient 

D = N dr dhiV^^'^l -^ 0,t)V''^'-^\l ^ 0,t - t)P, 

and a drift term 

-{Vr) drift = PeqlD 



Vr 



.dipeq 
) drift — h'eq'JJ^—T: • 



Explicitly, the diffusion coefficient has the form (see Appendix C^ and Ref. 



D = -NY^^rr'^^NPeJr). 
8 ' S r ^^ ^ 



;iii) 



:ii2) 



(113) 



(114) 



(115) 



;ii6) 



(117) 



(118) 
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3. The general case 

We can show in the general case that the relaxation of the test vortex is described by 
the Fokker-Planck equation 

^ + (V),,VP = V(D(VP + (3,,^PViJe,)) (119) 

with a diffusion coefficient 

^ = |^lniV(u;)e„ (120) 

where |S(r)| = 2J—Det(T,) is the local shear of the flow and and Det{T.) the determinant 
of the stress tensor S^'^ (see section |VC|) . In the regions where the shear cancels out, our 
approximations clearly break up. In particular, we cannot calculate the Kubo integral by 
assuming that the vortices follow the streamlines of the equilibrium flow. This is because, 
for a local solid rotation, the vortices always remain at the same relative distance and the 
correlation time is infinite. In that case it is necessary to take into account the dispersion 
of the vortices. An alternative derivation of the diffusion coefficient can be obtained by 
analyzing the statistics of velocity fluctuations created by a random distribution of point 
vortices |jl3[. When the differential rotation of the vortices is neglected (which corresponds to 



the opposite limit of that leading to equation ( |120[ )), we obtain (see Ref. |T^ and Appendix 

D^jVh^N (121) 

Clearly, a more complete study should take into account simultaneously the effect of the 
shear and the dispersion of the vortices to match the two formulae ( |120| ) and (|121|) . 



For t -^ +00, the distribution function P{r,t) of the test vortex converges towards the 



Maxwell-Boltzmann statistics (|26D . The time of relaxation corresponds typically to the 
time needed by the test vortex to diffuse over a distance R, the system size. Therefore 
treiax ~ R^ / D whcrc D r^ J lu N IS the order of magnitude of the diffusion coefficient given 
by equation ( [12CI| ). Using F = A^7 and introducing the dynamical time t^ ~ {uj)~^ ~ R^/^, 
we obtain the estimate: 

treiax ~ -, TT^-D- (122) 

m A 
Since the statistical description is expected to yield relevant results for large A^, we conclude 
that the relaxation of point vortices towards the Boltzmann distribution is a very slow 
process. It is plausible that a more violent relaxation be at work in the system. This 
problem is discussed more specifically in the conclusion. 

In the previous calculations, we have assumed that the distribution of the field vortices is 
given by the Maxwell-Boltzmann statistics ( PUj) which corresponds to statistical equilibrium. 
In the case of an arbitrary background distribution Peg, the expression of the drift is 

{V)dr^ft = DVlnPeg, (123) 

where D is still given by equation (|12UD. Since D > 0, the drift is always directed along the 



density gradient. The estimate of the time of relaxation is not changed in this more general 
situation. 
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V. A GENERALIZED KINETIC EQUATION 

In the previous section, we have described the relaxation of a test vortex in a "bath" 
of field vortices at statistical equilibrium. We would like now to relax this thermal bath 
approximation in order to obtain a generalized kinetic equation describing the evolution of 
the whole system. 

A. The factorization hypothesis 

If the vortices are initially decorrelated then, for sufficiently short times, they will remain 
decorrelated. This means that the (A^ + l)-particle distribution function can be factorized 
in a product of (A^ + 1) one-particle distribution functions: 

N 

fi{r,T,,...,T^,t) = l[P{rk,t). (124) 

fc=0 



If we integrate the Liouville equation ( |59D on the positions of the A^ vortices 1, ..., A^ and 
use the factorization ( |124|) , we directly obtain 

^ + (V)VP = 0. (125) 

Therefore, for sufficiently short times, the average vorticity (uj) satisfies the 2D Euler equa- 
tion. However, at later times, the distribution function fi differs from the pure product (124) 



and the Euler equation does not provide a good approximation anymore. In section [IV A| we 



have determined an exact equation ( p3D satisfied by the one-particle distribution function at 
any time. This equation is not closed, however, since it involves the A^-vortex distribution 
function fisys- We shall close the system by assuming that fisys can still be approximated by 
a product of N one-particle distribution functions: 

N 

/i^y,(ri, ..., Tat, t) ~ n ^(^k, t) (126) 

fc=i 

but, contrary to section [iV D| , the one-particle distribution function P(rfc,t) is not ascribed 
to the equilibrium value Peg. Physically, the decomposition ( |62D ( |126| ) assumes that the 
correlations that develop between the vortices (term /i/) are due uniquely to the polarization 
cloud imposed by each individual vortex. Without this polarization, the vortices would be 
uncorrelated (term figys)- In particular, this decomposition does not take into account three- 
body encounters which can play a crucial role in the dynamics of vortices (in particular when 
the system is neutral and homogeneous, see Ref. p7| , ^ ). These high order correlations 
develop on longer time scales and may be neglected in a first approach. 

The approximation (|126D introduced in equation (^) leads to a generalized kinetic equa- 



tion 

dP „ ., dP d 
'di 



N „t- N N 



dr drt^J ^Vi Jo ir[p^ 



8 8 \ ^ 

x^(t,t-r)(V'^(j-0)— + V'^(0-j)— jP(r,t-r)nP(r„t-r). (127) 
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Repeating the steps leading from equation (^) to equation ([99|), we obtain 

r f)P f)P ^ 

x{v^(l -. 0)Pi— + V^(0 - 1)P^}^_^, (128) 

where P = P(r, t) and Pi = P(ri, t). We also recall that between t and t — r the trajectories 
of the particles are determined from the smooth velocity field created by the vorticity distri- 
bution (uj) = N'yP{T,t). This non Markovian integrodifferential equation is similar to the 
equation obtained by Chavanis (2000) ||26| in two-dimensional turbulence using a quasilinear 
theory of the Euler-Poisson system. It is shown in Appendix that this equation rigorously 
conserves angular momentum in a circular domain and linear impulse in a channel (or in 
an infinite domain). However, under this form, it is not possible to prove the conservation 



of energy and the H-theorem. In our previous investigation of the problem [^, we have 
considered a somewhat crude approximation which amounts to neglecting memory effects 
in equation ( |128|) . We beleived that this approximation would not introduce any significant 
error but, as we shall see, we were wrong: such approximation breaks the conservation of 
energy. 

If we assume that the decorrelation time r is short (which does not need to be the case) 
and implement a strong Markov approximation, we obtain 

In the case of an infinite domain V(0 ^ 1) = — V(0 -^ 1) and we have the further simplifi- 
cation 

dP , ,^,,dP N-f\ d f. r.fuu^^^fr. 9P ^dPi 



where 



IF + <^> 97 = ^^*? / '"■'^^'""® l^'*^ - ^*f )• 'i^"' 



K'^^m = -^ = ^^^ (131) 

and ^ = ri — r (we also recall that ^_|_ is the vector ^ rotated by +f )• To arrive at equation 
( |131| ) we have explicitly used the form of the Kernel (|^ and to get the second equality 
we have used the fact that we are in two dimensions. Note that the symmetrical form 
of equation ( |130| ) is reminiscent of the Landau equation introduced in plasma physics and 
in stellar dynamics |0]. In this analogy, the position r of the vortices plays the role of the 
velocity v of the electric charges or stars. Therefore, we can directly infer the conservation of 
linear impulse P_l = J urd'^r and angular momentum L = J ur'^d'^r which play respectively 
the role of impulse P = J fvd^v and kinetic energy K = ^ f^d^w in plasma physics. 
In addition, we can prove a H-theorem for the Boltzmann entropy (P^ ) exactly like for the 
Landau equation. Finally, we can show that the solutions of equation ( |130| ) converge towards 
the Gaussian vortex (the equivalent of the Maxwellian distribution in plasma physics with 
r in place of v): 
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■iQ7(r-ro)^ 



(132) 



which is the maximum entropy state at fixed circulation, angular momentum and impulse. It 
is in general different from the Boltzmann distribution (p6|) with the relative streamfunction 
ip' = ip + ^r^ — U_L.r. This clearly indicates that equation ( |129| ) does not conserve the 
meanfield energy (0). It may happen, however, that the energy is approximately conserved 
if we start from an initial condition with a value of energy Eq corresponding to /? — >■ 0, 
Q -^ +00 and a = (3VL/2 finite at equilibrium. In that case, equation ( |132| ) is the maximum 
entropy state at fixed F, L and Eq (the conservation of impulse can be satisfied trivially 
by taking the center of vorticity as origin of our system of coordinates). However, if E 
differs from Eq by a large amount, the kinetic equation (|129|) will not correctly describe the 
evolution of the system for late times. 

Now, if we account properly for memory effects in equation ( |12^ ), we can obtain a more 
general equation which guaranties in addition the conservation of energy and is therefore 
more satisfactory. In the case of an axisymmetrical flow, it is possible to calculate the 
memory function appearing in equation (|128|) explicitly if we assume that the correlation 



time is smaller (but not necessarily much smaller) than the typical time on which the average 
vorticity changes appreciably. In this approximation, the point vortices follow, between t 
and t — r, circular trajectories with angular velocity r2(r, t) = {Vg){r, t)/r and equation ( |128| ) 
simplifles in (see Appendix ^: 



dP _ N-f^ d 
dt Ar dr 
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(133) 



where Vt = Q{r,t), Qi = Q{ri,t) and r> (resp. r<) is the biggest (resp. smallest) of r and 
ri. The angular velocity is related to the vorticity by 

1 d 



(uj) 



r dr 



{nr' 



(134) 



Similarly, in the case of a unidirectional flow, we can assume that, between t and t — t, the 
point vortices follow linear trajectories with velocity (V) = {V){y,t)ex. This leads to the 
kinetic equation (see Appendix ^: 



dP 
'dt 



~^dy 



+ 00 
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(135) 



where V = {V){y,t) and Vi = {V){yi,t). The function -E'i(x) is the exponential integral and 
L an upper cut-off necessary in that case (see Appendix |^). We also recall that the average 
velocity is related to the vorticity by 



(.) - -|w. 



(136) 



We can remarkably propose an approximation of the general kinetic equation ( |128| ) which 
encompasses both the axisymmetric form (|133|) and the unidirectional form (|135|) . Memory 
effects are not neglected, unlike in Eq. ( |130| ), but they are simplifled in a way which preserves 
all the conservation laws of the system (as discussed below). We propose the generalized 
kinetic equation: 
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with 

K^\i) = -^ = ^^^ (138) 

and ^ = ri — r, V = (V)(ri,t) — (V)(r, t). For specific applications, it may be necessary to 
introduce a shielding of the form ([E39|) in the interaction between vortices. This shielding 
arises naturally in geophysics in the "quasigeostrophic approximation". In that case the 
tensor K^" is replaced by (1/L)ii'i(^/L)^^^^/^ where K^ is the modified Bessel function of 
first order and L is an upper-cut off (called the Rossby radius in geophysics) which plays 
the same role as the Debye length in plasma physics. 



Equation ( 137 ) is not exact, in a strict sense, yet it satisfies all the conservations laws 



of the vortex system as well as a H-theorem (section |V lj| ). This is very gratifying and 



this can have important practical applications. It is remarkable that we can write down 
an approximate kinetic equation in the general case without being required to specify the 



trajectory of the point vortices between t and t — r. In fact, to arrive at equation ( |137| ) we 
have made implicitly two approximations: (i) we have assumed that the vorticity field does 
not change dramatically when we follow the vortices in their motion between t and t — r. 
(ii) After the time integration has been effected, the non universal function V^{\ -^ 0,t — T) 
gives rise to a logarithmic term which has been replaced by 1 in the subsequent calculations. 
This term produces a sub-logarithmic correction that is flow-dependant and that has been 
neglected. It is on account of this weak dependance that a general kinetic equation can be 
obtained. For axisymmetrical or unidirectionnal flows, equation (|137|) reduces to equations 



( |133|) and ( [L35|) with ln(|ri — r|/(ri -|- r)) and Ko{\y — yi\/L) instead of ln(l — (r</r>)^) 



and Ei{2\y — yi\/L). This is not a too severe discrepency (these functions have a similar 
logarithmic behaviour) so our approximations are reasonable. 

From equations ( |133| ) ( |135| ) and ( |137| ), it is clear that the relaxation towards equilibrium 
is due to a phenomenon of resonance. Only the points ri satisfying the condition^.v = with 
ri 7^ r will contribute to the diffusion current in r. In the axisymmetrical case, this condition 
of resonance reduces to fi{ri) = fi{r) and in the unidirectional case to V{yi) = V{y). These 
conditions of resonance had never been noticed previously. Further work on the subject 
will have to make these criteria more precise by computing explicitly "resonance lines" in 
two-dimensional real flows. 



B. Conservation laws and H-theorem 

We now derive the conservation laws and the H-theorem satisfied by equation ( |137D 



First of all, the conservation of the circulation is straightforward since equation (|137|) can 



be written in the form of a continuity equation ( pi|) . To prove the conservation of angular 
momentum, we start from equation ( p4| ), substitute for ( p.37| ), permut the dummy variables 
r and ri and add the resulting expressions. This yields 

L = ^Jd'-dh,K^^em-){pi^ - p^). (139) 
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But, from equation (|138|) , we immediately verify that 

K^^'i^ = 0. (140) 

This proves the conservation of angular momentum. We can prove the conservation of 
linear impulse in a similar manner. Starting from equation (|D16|) , substituting for (|137D , 
permuting the dummy variables r and ri and adding the two resulting expressions yields 
P = 0. For the conservation of energy, we start from equation (|5^) and follow the by-now 
familiar procedure. This yields 

^"^/'''■■'''■■■*''""Wev)(F,^-F||). (141) 

Considering the form of the tensor ( p.38|) , we have 






^"X = 7^(ev). (142) 



When substituted in equation ( [L41|) , we see that the occurence of the delta function in 
the kinetic equation implies E = 0. Finally, for the rate of entropy production we have, 
according to equations ( ^31) and (|137|) : 

8 J PPi dr^ ^^ \ dr" dr'( J ^ ^ 

Permutting the dummy variables r and ri and adding the resulting expression to equation 
( |143|) , we get 

16 J PPi ^^ \ drt^ dr^J V dr'' dr'( J ^ ' 

Now, for any vector, A^^K^'^A" = (A|'_l)^/^^ > 0. This proves a H-theorem (S > 0) for 
our kinetic equation ( |1371) . It should be emphasized that the conservation laws and the H- 
theorem result from the symmetry of the kinetic equation (and the condition of resonance) 
and not from formal Lagrange multipliers like in the thermodynamical approach of section 
[111 C| . In addition the H-theorem is proved by our approach instead of beeing postulated. 
This is more satisfying on physical grounds. 

It remains for one to show that the Boltzmann distribution 

P = ^e-/37(^+|r2-Ux.r) (145) 

is a stationary solution of equation (|137|) . Noting that 

|^ = -/57(^ + ^-^-f/l)i^, (146) 

we have successively 

/ dP dP \ f^ 

K'\PiQ^ - P^j = PlPPiK^^ivl + nn = PlPP.fiC^), (147) 
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where we have used equations (|140|) and ( |142| ). When substituted in equation (|137|) , we 



find that the r.h.s. cancels out due to the deha function. The advective term is also zero 
since P = f{ip')- Therefore, the distribution ( |145| ) is a stationary solution of equation ( |137D . 
Note, however, that this is not the only solution. Any stationary solution satisfying ^.v 7^ 



for any couple of points r, ri (with r 7^ ri) is a solution of equation ( |137| ). Physically, this 
implies that the system needs sufficiently strong resonances to relax towards the maximum 
entropy state. If this is not realized it will be frozen in a sort of "metastable" equilibrium. 
This may explain why the maximum entropy state is not always reached in two-dimensional 
turbulence. For example, for a unidirectional fiow with (u) positive or negative everywhere, 
the velocity field is monotonous (see equation ( |136| )) and the condition of resonance cannot 
be satisfied. The evolution of the system will require non trivial correlations between point 
vortices that are not taken into account in the present approch. One would need to replace 
the factorization hypothesis ( |126| ) by a product of two-point or three-point correlations 
functions. However, it is plausible that these correlations develop on a very long time scale 
so it remains a matter of debate to decide whether they really are relevant for the dynamics 
or not. 

In the context of 2D turbulence described by the Euler-Poisson system, the quasilinear 
theory developed by Chavanis (2000) yields instead of equation ( [1371 ): 

duJ _ e^ (9 /■ ,0 ,_„._.^ . r ,, ,. duJ , .^uJ'^ 

i4s 



^ + "^^ = ¥^ / ^'r'i^-5(ev){^Vo - -0^ - ^(-0 - ^) 



^r'^'}' 



where K^^'^ is defined by ( |138|) and uJ = cJ(r, t), uJ' = cJ(r',t). In addition to the previ- 
ous conservation laws, this equation guaranties that the coarse-grained vorticity uJ remains 
bounded by the maximum value of the initial distribution, i.e. cJ < ctq. This equation 
satisfies a H-theorem for the Fermi-Dirac entropy introduced by Miller-Robert-Sommeria at 



equilibrium p8|j29|| . Our approach provides therefore another way of justifying their results 
from a dynamical point of view. Equation (|148|) is written for a single level of vorticity ctq, 
but it is possible to extend the quasilinear theory to an arbitrary distribution of vorticity 
levels (in preparation). Our equations should provide therefore an interesting and useful 
parametrization of the 2D Euler equation. It should be recalled in that respect that the 
usual turbulent diffusion z/AU introduced ad hoc in the r.h.s. of the 2D Euler equation in 
order to smooth out the small scales and prevent numerical instabilities breaks the conser- 
vation laws of the inviscid dynamics. This is not the case for our equation ( |148|) : not only 
it smoothes out the unresolved scales (as exemplified by the existence of a H-theorem) but 
it satisfies all the conservation laws of the inviscid dynamics and respects the invariance 
properties of the Euler equation (invariance by translation and rotation of the coordinates, 
Galilean invariance and invariance by rotation of the referential). In addition there is no 
free parameter in our theory except the coarse-graining mesh e (or resolution scale) which 
depends on the situation contemplated. Different attempts had been made previously to 
obtain an equation satisfying all these requirements, but only partial results were obtained 

il,liiij2i. 
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C. Connexion ^vith the Fokker-Planck equation 

Equation (|137|) can be considered as our final result but we wish to show that a direct 
connection with the Fokker-Planck equation of section ^ can be found. Introducing a 
diffusion tensor 



jjf^u ^ ^_ j dhiK^''5{^.v)Pi (149) 

8 J 



and a drift term 



^. = _N^' f.2,- r.u.......dP, 



J dh.K^'^Sitv)^, (150) 



equation (|137|) can be rewritten in the more illuminating form: 

(151) 



dP d 

+ (V)VP 



dt drf^ 



dP 



similar to the general Fokker-Planck equation (^Ol). Note, however, that equation ( |151| ) is an 
integrodifferential equation since the density probability P(r, t) in r at time t depends on the 
value of the whole distribution of probability P(ri, t) at the same time by an integration over 
ri. By contrast, the Fokker-Planck equation (^Tj) is a differential equation. The usual way 
to transform an integrodifferential equation into a differential equation is to make a guess 
for the function -P(ri) appearing under the integral sign and refine the guess by successive 
iterations. In practice we simply make one sensible guess. Therefore, if we are close to 
equilibrium, it seems natural to replace the function Pi appearing in the integrals by the 
Boltzmann distribution 

P(ri) = Ae-^^^'^'^\ (152) 



This corresponds to the "thermal bath approximation" of section P^VD| : the vortices have 
not yet relaxed completely, but when we focus on the relaxation of a given point vortex 
(described by P) we can consider, in a first approximation, that the rest of the system 
(described by Pi) is at equilibrium. Within this approximation, the diffusion coefficient and 
the drift simplify in: 

di)' 

D^" = ^Pir, t) J K^''6{t^r)d% (154) 

where we have made the local approximation. If we assume that the correlation time is 
short, i.e. if we replace C^5(^.v) by t/tt^ (compare equations ( |1371 ) and (|130|) ), we obtain 

ri = p-fDVip', (155) 

D = ^\nN{uj). (156) 
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In that case, equation (|151|) reduces to the Fokker-Planck equation found in section fV : 



dP 

-— + {V)VP = V{D{VP + p-fPV^')). (157) 

ot 

If instead of the Boltzmann distribution (|1 52|) we use the Gaussian distribution (|132D , we 
get 

dP 

— + (V)VP = V(D(VP + a7Pr)). (158) 

This equation is closely related to the Kramers- Chandrasekhar equation (|^) since the drift 
and the friction are linear in r and v respectively. 

The diffusion coefficient (|156|) was previously obtained by [^^,^ using phenomeno- 



logical arguments. However, in these studies the correlation time r was not specified. A 
first determination of r was obtained in Ref. 0] but it was restricted to axisymmetrical or 
unidirectional flows. Using equation (|154|) , we can determine the expression of the diffusion 
coefficient and the correlation time in the general case. Expanding the velocity difference 
V = (Vi) — (V) in a Taylor series in ^ = ri — r, we obtain to ffist order in the expansion 

e.v = S^^e'^C, (159) 

where 

E- = if^ + ^) (160) 

is the stress tensor. It satisfies the property of symmetry S^'^ = S"^^. Since the flow 
is divergenceless, we also have T,^^ + S^^ = 0. This suggest to introduce the notations 
a = E^^ = — S^^ and b = T/^'^ = S^"^. In terms of the stress tensor ( |160| ), the diffusion tensor 
( |154|) can be rewritten 

O J (^ 

This integral can be performed easily by working in a basis where the tensor S'^'^ is anti- 
diagonal. To that purpose, we seek a tensor T such that S = TH'T where S' is anti- 
diagonal and, by definition, T^'^ = T^^. We also impose that T is unitary so that T = T~^. 
Then, if we denote by (Af-*^^, M-*^^, M^^, M^^) the components of a 2 x 2 matrix, we find 
S' = (0, h\ h\ 0) and T = (a, /5, -/3, a) with a^ + /^^ = 1 and V = b/{a^ - l3^) = -a/2a(3. 
From the above results, it is also clear that 6'^ = — Det(S) = (a^ + 6^), where Det(S) 
stands for the determinant of the matrix S. Now, introducing a new system of coordinates 
such that ^ ^ = T'^'^^'^, or alternatively ^^ = T"'^^ '', we easily check that the Jacobian of the 
transformation^ -^ ^, i.e. the determinant of T, is equal to one. Under these circumstances, 
the diffusion tensor (|161|) can be written D = TD'T with 

D'"" = -^pj - — ^^^K\n^mQ de,d^'„ (162) 



where we have set |S(r)| = 2J —Det{Y?). Physically, this quantity represents the local shear 
of the flow. The components of the tensor (|162|) can now be determined easily. First of all, 
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D"' = -^T^^P / -Jh^mQ di[di', (163) 



" |S(r)| J e? + e? 
Setting i'i= i cos ^ and ^'2 = i sin ^ where ^ = ^' = |ri — r|, we get 

TV^^ 1 /■+00 /•27r 

^^"= : 1^, ,, P idi d9sm'9S{ecosesm9), (164) 

8 |z^(r)| JO Jo 

or, equivalent ly, 

^,11 ^ ^J ^ p f+'-rf^ r desmOdicose). (165) 



4 |S(r)| Jo e ^0 

As explained in Appendix ^ we regularize the logarithmic divergence by introducing ap- 
propriate cut-offs at small and large scales. With the change of variables t = cos 9, we 
obtain 



By the same arguments, we find that D'"^^ = D'^^ = D. On the other hands, it follows for 
reasons of antisymmetry by the transformation ^[ —>■ —^[ that D'^^ = D'^^ = 0. Therefore 
D'fj-'^ = ]J5^'^ is diagonal in the basis where S is anti-diagonal. This remains true in any 
basis since D^''' = T^t^D'^'^T'"' = D{fTy = Dd^"" . Therefore, close to equilibrium, the 
diffusion is isotropic and the general expression of the diffusion coefficient is 

Comparing with equation ( p.56| ) we find that the correlation time is given by 



The quantity |S(r)| = 2J —Det(Yi) plays a fundamental role in the theory. Clearly, this ex- 
pression is invariant by a change of referential. For a unidirectional flow or an axisymmetrical 
flow, we recover the results of section |1V \^ and of Ref . |p . 

VI. CONCLUSION 

In this paper, we have provided a systematic derivation of the kinetic equations of point 
vortices, applying for the first time the powerful projection operator technics to this problem. 
We have described how a cloud of point vortices relaxes toward the mean field statistical 
equilibrium, leading to a clustering into large coherent vortices. In the first part of the paper, 
we have focused on the relaxation of a "test" vortex in a cloud of background vortices at 
statistical equilibrium. The cloud of "field" vortices plays the role of a thermal bath like in 
other problems of statistical physics. We have shown that the test vortex undergoes a usual 
diffusion effect due to random fluctuations and that it also experiences a systematic drift. 

31 



This drift, due to a polarization of the background vortices by the test vortex, balances the 
effect of diffusion at equilibrium, providing a dynamical explanation for the persistence of 
clustering. The drift was previously derived with a linear response theory (Chavanis, 1998), 
but the diffusion was heuristically introduced by adding a white noise effect. The present 
derivation systematically derives the two effects, diffusion and drift, from the same formal- 
ism, and is therefore more satisfactory. The diffusion derived here turns out to be influenced 
by long time correlations, so it is more complex than the usual white noise effect. When 
memory effects are ignored, we obtain a Fokker-Planck equation for the evolution of the 
one-particle distribution function. This Fokker-Planck equation can also be derived from a 
phenomenological Maximum Entropy Production Principle ||2^. This shows that the struc- 



ture of this equation is influenced more by thermodynamics (the first and second principles) 
than by the precise microscopic model. However, our systematic procedure starting directly 
from the Liouville equation provides a justification for this thermodynamical approach and 
specifies its range a validity. It also allows to determine explicitly the value of the diffusion 
coefficient which was left unspecified by the Maximum Entropy Production Principle. All 
these results could be tested numerically by introducing a test vortex in "sea" of vortices at 
statistical equilibrium and by solving the Kirchhoff-Hamilton equations of motion. 

In the second part of the paper, we have attempted to describe the evolution of the whole 
system of vortices far from equilibrium. We have obtained a new kinetic equation (|128|) 
which incorporates a delocalization in space and time. This is therefore a non Markovian 
integrodifferential equation. A similar equation also occurs in the quasilinear theory of the 
2D Euler equation (Chavanis, 2000). Within some approximation, it is possible to carry out 
the time integration explicitly and this yields a simpler equation ( p.37| ) which only conserves 
a delocalization in space. This equation respects all the conservation laws of the point vortex 
system and satisfies a H-theorem. The relaxation is due to a condition of resonance between 
distant vortices. If the system is sufficiently "resonant" , it will reach a maximum entropy 
state described by the Boltzmann distribution. However, if there are not enough resonances 
the evolution may stop on metastable state. Only non trivial correlations between vortices 
(for example three body collisions) can unfreeze the system and induce further evolution. 
These correlations are not taken into account in the present analysis although the projection 
operator formalism might still apply. It would be necessary to modify the factorization 
hypothesis ( |126| ) so as to account for two-body or three-body correlation functions. It is 
possible, however, that these correlations develop on a much longer time scale so it is not 
yet clear whether they are physically relevant. In any case, the approximations made in the 
present paper are a first step towards a rational kinetic theory of point vortices. 

It is also possible that a system of point vortices undergoes a form of "violent relaxation" 
in its early stage. For short time scales, the correlations between point vortices have not yet 
developed and the average vorticity is solution of the 2D Euler equation. When the initial 
condition is far from equilibrium, it is well-known that the 2D Euler equation develops a 
complicated mixing process leading to the formation of an organized state (on a coarse- 
grained scale). This relaxation is quite rapid, of the order of the dynamical time t£,, and 
the resulting equilibrium state is predicted to be a complicated superposition of Fermi- 
Dirac distributions respecting all the constraints of the Euler equation |2^J2^. On longer 



time scales, of the order of ^^to, the correlations between point vortices develop and 
the system undergoes another form of relaxation, much slower. This relaxation is towards 
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the Boltzinann distribution derived by [^ Joyce & Montgomery and [Q Lundgren & Pointin 
which is the true equihbrium state for a system of point vortices. The first type of relaxation, 



by vorticity mixing, was described in Ref. p6| using a quasihnear theory of the 2D Euler 



equation. The exclusion principle leading to the Fermi-Dirac statistics was explicitly shown 
as well as a if theorem. The second type of relaxation, due to discrete interactions between 
point vortices, was the object of the present paper. 

It is noteworthy that a similar distinction occurs in the context of stellar systems 0. 
Indeed, the relaxation of stars is a two stage process. For short time scales ~ to, the 
encounters between stars can be neglected and the distribution function is solution of the 
Vlasov equation (analogous to the 2D Euler equation). If the system is initially far from 
mechanical equilibrium, it will experience a "violent relaxation" towards a virialized state. 



This equilibrium is predicted to be a superposition of Fermi-Dirac statistics [pT| , p2| , like for 



the 2D Euler equation. Then, on a longer time scale, of the order of i^^d, the encounters 
between stars cannot be ignored anymore and will deviate the stars from their unperturbed 
trajectories. This coUisional relaxation is usually described by a Landau or Fokker-Planck 
equation that converges towards the Maxwell-Boltzmann distribution at equilibrium. On 
even longer time scales, three-body encounters leading to the formation of binaries induce a 
slow evolution of the system |^3[. Three body encounters (involving vortices of different sign) 



are also relevant in two dimensional turbulence and lead to the formation of vortex "pairs" 



2^ . The analogy with "binary stars" is interesting to note. It would be important to test 
these ideas with numerical simulations of point vortices or stars. The situation is difficult 
in the stellar context because a maximum entropy state does not always exist. Indeed, the 



system can collapse and overheat: this is the so-called "gravothermal catastrophe" |33]. This 



problem does not occur for point vortices and it should be possible to evidence the two (or 
more) successive equilibria more properly. An advantage of point vortices with respect to 
stars is the lower dimensionality of space [D = 2 instead oi D = 3, or D = 6 in phase space) 
that should make numerical simulations easier. 
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APPENDIX A: THE STATISTICS OF VELOCITY FLUCTUATIONS IN AN 

INHOMOGENEOUS MEDIUM 

In this Appendix, we study the statistics of velocity fluctuations produced by an inho- 
mogeneous distribution of point vortices. This study extends the calculations of ||35| , |36|jr3[ 



for a uniform medium and provides a simple framework to understand the logarithmic diver- 
gence of the diffusion coefficient. Let us consider a collection of A^ point vortices randomly 
distributed in a disk of radius R with an average density nij). The velocity V occuring at 
a given location r of the flow is the sum of the velocities $i {i = 1, ...,N) produced by the 

N vortices: 

V = f:*.. (Ai) 



*. = -f f^. (A2) 

2ti |rj — rp 

Following a procedure similar to that adopted in Ref. [1^, the velocity distribution can be 
expressed as 

Wn{Y) = ^lAr,ip)e-^d'p (A3) 

with 

A^ip) = I^Je^P{n)dh,y. (A4) 

Here, -P(ri) denotes the probability of occurence of a point vortex in ri and by definition 



$= — -^ ^. (A5) 

27r ri -r 2 ^ ^ 



Introducing explicitly the vortex density n{ri) = NP{ri), we have 

A^ip) = (^^ J n{r,)e^dh,y . (A6) 

Since 

f n{ri)dhi = N, (A7) 

we can rewrite our expression for An{p) in the form 

AM = (l - ^/n(rO(l - ef^)dh,y. (A8) 

In the limit of large A^, R with fixed n(r), we can approximate the foregoing expression by 

A{p) = e-^^^ (A9) 
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with 

C{p) = fn{r^){l - e^)t/2ri. (AlO) 

Separating the real and imaginary parts of C{p), we obtain 

C{p)=C,{p)-zC2{p) 
= / n(ri)(l - cos(^))(i^ri - i / n(ri) sm{f^)(fri. (AH) 

In the first integral, we find it convenient to introduce the relative separation ^ = ri — r in 
terms of which 

Ci(p) =|n(r+e)(l-cos(^))rf2e (A12) 



with 






* = -T^tI- (A13) 



In Ref. |TB|, it was found that an important contribution to the velocity fluctuations comes 



from the nearest neigbour. This justifies to make the local approximation n(r +^) ~ n(r) 
in equation ([A12D . In this approximation 



Ci(p) = n(r) [{1 - cos(/^))d% (A14) 



This quantity is closely related to the function C{p) evaluated in Ref. |T^ for a distribution 
of vortices with uniform density n. The only difference is the presence of the local density 
n(r) in place of n. Therefore, we can infer directly that 



C.W = ^^1„^— )/. (A15) 



In fact, the local approximation is only marginally valid because, as discussed in Ref. JTS 
the contribution of the nearest neighbor is precisely of the same order of magnitude as the 
contribution of the rest of the system. This results in a logarithmic divergence in ( |A15| ) due 
to the weak collective behaviour of the system. 



In the second integral appearing in ( [All| ), the contribution from proximate vortices 



vanishes by symmetry. As a result, the integral is dominated by large values of |ri| or, 
equivalently, by small values of |$|. We can therefore make the approximation sin(/5^) ~/5^ 
and write 

C^ip) = p j n{vr)^d\r. (A16) 

In the integral, we recognize the mean-field velocity created in r by the average distribution 
of vortices: 

{V){r)= fn{ri)^dhi. (A17) 
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Hence, 

C,{p)=p{Y){v). (A18) 

Substituting the explicit expression for C{p) in equation ( |A9|) , we obtain 

A(p)=e-^^'^(^)''^+W(r). (A19) 

Therefore, the velocity distribution in r can be written quite generally 

Wiy) = ^J e-^ W^)P^e-^(v-(v)(r))^2^_ (A20) 



This is the same distribution as for a uniform distribution of vortices except that the constant 
density n has been replaced by the local density n(r) and that the distribution is for the 



fluctuating velocity V = V — (V)(r). Repeating the calculations of Ref. ||T3[, the velocity 
p.d.f is explicitly given by: 

4 47r ■,;2 

^^^^ = n(r)yiniv """^'^^^ (V < V,„,(iV)) (A21) 

^(^)-J^ (V>Ve„,(Ar)) (A22) 

where 

VcritiN) ~ (^!^^\nNY\n'/\\nN). (A23) 

This distribution lies at the frontier between Gaussian and Levy laws: the core of the distri- 
bution is Gaussian while the tail decays algebraically like for Levy laws. This is because the 
variance of the individual velocities ( [A2|) diverges logarithmically so the central limit theo- 
rem is only marginally applicable. For that reason we have proposed to call this distribution 
the "marginal Gaussian distribution" |]14[. In the strict mathematical limit A^ —>■ +00, the 
transition between the two regimes is rejected to infinity and the velocity p.d.f. is purely 
Gaussian [Q. However, the convergence towards this Gaussian distribution is so slow that 



in practical applications it is never reached: the algebraic tail always remain pT |. 



According to equations ( |A1| ) and ( |A2| ) , the variance of the velocity can be written 



where we have made the local approximation in the second equality. This quantity diverges 
logarithmically at both small and large vortex separations. The divergence at small sepa- 
rations is a failure of our model which ignores correlations between vortices. In fact, when 
two vortices approach each other they can form a pair, as discussed in Ref. [21], and our 
mean field theory clearly breaks down. We shall account heuristically for this failure by in- 
troducing a cut-off at some minimum distance dpair ~ (vrnln A^)~^/^ [|l^. The divergence at 



large separations is due to the unshielded nature of the interaction potential. It is therefore 
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natural to cut the integral at i?, the typical size of the system. With this regularization, we 
obtain 



(VO 



2\ ^(r)7% fR\ «(r)7^ 



2n ^"UJ 



An 



\nN. 



(A25) 



Since the divergence is weak (logarithmic), the result does not depend crucially on the precise 
value of the cutoffs. The same expression for the variance can also be obtained from the 
formula 



(V^) = f W{V)V^2TTVdV 



(A26) 



if we introduce a cut-off at large velocities and use equations ( |A21| ) ( |A22|) ([A23| ) [to sufficient 
accuracy, we just need considering equation ( [A21| )]. 

The diffusion coefficient can be expressed in terms of the variance ( |A25|) and the typical 
correlation time r by [ll3| : 



D = -^r{V'). 



Using equation ( A25 ) and the relation (uj) = wy, we obtain 

'-fT 



D 



IGvr 



IniV(cj) 



(A27) 



(A28) 



in agreement with our result ( |15tj| ). The correlation time is more difficult to evaluate. If 
we ignore the differential rotation of the vortices, the calculations of Ref. |]13| are directly 
applicable and lead to the expression 



1 



r 



(a;)(r)VIKiV' 



(A29) 



'n r 



This is the typical time needed by a vortex to cross the interparticle distance rf ~ l/wr 

(V^). With this approximation, the diffusion coefficient is given by 



with the velocity 

equation ( |121| ). On the other hand, if we consider that the vortices are transported by the 
equilibrium flow and evaluate the diffusion coefficient with the Kubo formula (see Appendix 
P), we find that the correlation time is related to the local shear by equation (|168|), i.e.: 



27r 



\nr) 



(A30) 



Physically, it corresponds to the time needed by two vortices with relative velocity Sd to be 
stretched by the shear on a distance ~ d. This approximation breaks up, however, when the 
shear is weak. In that case r is given by equation ( |A29D obtained when only the dispersion of 
the vortices is considered. Clearly, a general formula should take into account simultaneously 
the effect of the shear and the dispersion of the particles. 
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APPENDIX B: THE KUBO FORMULA 

Let us consider the diffusion of a test vortex in a "sea" of field vortices described by the 
equilibrium distribution 

TV 

/ieg(ri, ..., rjv) = n Peq{ri). (Bl) 

The general form of the diffusion coefficient, as defined by equation (ffDI), writes 



Now, the net displacement of the test vortex produced by the fluctuations of the velocity 
between t and t + At is given by 

ft+At 

Ar = / Vit')dt'. (B3) 



Substituting explicitly for Ar from equation ( p3D in equation (|B2|) , we have 



T fAt i-At 

^'^^2AtJo "^^'i dt"{V''it + t')V''it + t")),,, (B4) 

where {.)eq denotes the average with respect to the equilibrium distribution ( pi]) . Since the 
correlation function appearing in the integral only depends on the time difference \t" — 1'\, 
we also have 

L>'^-= / dt' dt"{V^{t + t')V''{t + t")),, (B5) 

At JO Jo 

or, alternatively, 

D^^ = — dt' / t/t"(V'^(t)V''(t + t" - t'))en. (B6) 

At JO JO 

With the change of variables t = t' — t" , we obtain successively 

1 /-At fAt 

If the correlation function (V'^(t) V(t — T))eq decays more rapidly than r^^, we can take the 
limit At -^ +CX) to finally obtain 



r+oo 

^M- = / (V^(t)V'^(t - r))e,rfr. (B8) 
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This is the Kubo formula for our problem. Remembering that V denotes the fluctuation of 
the total velocity: 

V(t) = V(t) - (V(t))e„ (B9) 

we find that 

{V^{t)V^{t - r))e, = i,V^{t)V\t - r))e, - {V^(t)),,iy\t - r)),,. (BIO) 

Now, the first quantity in bracket can be written explicitly 



N „ N 

= E E y V'i^ ^ 0, t)\/^(j ^ 0, t - r)/ie,({r J) n rf'rfc 

j=l jyj k=l 

N „ N 

+ E y ^^(^ ^ 0' ^)^'(^ ^ 0, t - r)/ie,({rfc}) n dhk 



i=l k=l 

= N{N - 1) {V^il -> 0, t)),g{V''{l ^ 0, t - r))e, 
+iVy'r'^(1^0,t)r"(1^0,t-r)Pe,(ri)rfV. (Bll) 



For large A^, we can make the approximation A^(A^ — 1) ^ A^^. Using ( |B1(J| ) and ( |B11| ) the 
correlation function can be put in the form 

{V>^{t)V'{t - r))e, = nJ \/'^(1 -^ 0, t)V''{l ^ 0, t - T)P,,ir,)dh,. (B12) 

Since the integral is dominated by interactions involving relatively close vortices, we can 
make the local approximation: 

(V'^(t)V^(t - r))e, = NJv^il -^ 0, 1)^(1 ^ 0, t - T)P,^{r)dh,. (B13) 

The expression for the diffusion coefficient then becomes 

D^"' = N dr ¥"{1^0, t)V%1^0,t- T)P,g{r)dhi. (B14) 

For sake of brevity, we shall denote the velocity correlation function by 

C^'^ir) = {V^{t)V'{t - r))e, = nJ V^{1 ^ 0, t)^^! ^ 0,t - r)P,g{r)dh,. (B15) 

Therefore, the Kubo formula takes the form 

£)M- = / C'"'{T)dT (B16) 

Jo 
More generally, we have 

/•At 

(Ar^Ar'') = 2 / C^^(r)(At - T)dT. (B17) 

Jo 

These quantities are calculated explicitly in Appendix P in the case of simple flows. 
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APPENDIX C: THE CALCULATION OF THE DIFFUSION COEFFICIENT 

In this Appendix, we calculate the Kubo integral using an approximation in which the 
point vortices follow the streamlines of the equilibrium flow. 

1. Unidirectional flow 

We shall first calculate the velocity correlation function ( |B15| ) and the diffusion coefficient 
( pi4|) in the case of an unidirectional equilibrium flow. The trajectory of a fluid particle 



advected by this flow is simply: 

y(t-r)=y{t), (CI) 

x{t-T) = x{t)-ly),,{y)r. (C2) 



According to section |I V E 1| , we are particularly interested by the yy component ( |103| ) of the 



velocity correlation function. Explicitly, it has the form: 

where we have used equation (|^). The second term involves the quantity 

(xi - x)(t - r) = xi - X + {{V),,iy^) - {V),,{y))T. (C4) 

In the local approximation, we can expand the velocity difference in a Taylor series in yi —y. 
To first order, we have 

{V)e,iyi) - {V)e,iy) ^ -ny){yi - y). (C5) 

where S(|/) is the local shear of the flow ( |1U5| ). Introducing the variables X = xi — x and 
F = l/i — y, we obtain: 

The integration over X can be performed easily since the integrand is just a rational function 
of polynomials. After straightforward calculations, we find: 



^p <y) \ P^ 



C{r) = -^PeM . , 1^...,. /„ — • (C7) 



The integral over Y diverges logarithmically for both small and large Y . The reason for 
this divergence has been given in Ref. |T^ and in Appendix 0. Introducing two cut-offs at 
scales d and i?, and noting that ln(i?/(i) ~ ilnX, we finally obtain equation ( |104D . For 



T -^ +00, the correlation function decreases like r ^. This is a slow decay but still the 
diffusion coefficient (|B16| ) converges. Using 
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(C8) 



and taking the limit t -^ +00, we find equation (|109|) . More generally, using ( P17|) and 
^81), we have 



{{^yf) 



N-f"^ InN 



PeM 



2vr \ny)\ 
For At — i> (ballistic motion): 



arctanQ|S(y)|At)At - -i-ln(l + h:{yf{^tf) 



■ (C9) 



N^^ 



{{Ayr) = -^\nNP,,{y){Aty = -{V'){Atr 



in 



1 

2' 



and for At -^ +00 (diffusive motion): 



((A,)^) = ^i^ 



Peg{y)At. 



(CIO) 



(Cll) 



2. Axisymmetrical flow 

In an axisymmetrical flow, the trajectory of a fluid particle takes the simple form: 

r{t-T)=r{t), (C12) 



9{t -t) = d{t) - ^^-h^ 



T. 



(C13) 



As indicated in section IV E 2 , we are particularly interested by the r{t)r{t — r) component 
( |112| ) of the correlation function. Let us introduce the separation 5r = ri — r between the 
field vortex 1 and the test vortex. In the local approximation, 5r can be considered as a 
small quantity. Therefore we can write: 



5r = rSOeg + 5re,. = Xeg + Fe^, 



(C14) 



dhi = cf{5r) = dXdY. 

With these notations, the correlation function ( |112D can be rewritten: 

N'y'^ r X X 

C{t) = -\Pe,{r) / dXdY———{t)———{t - r). 



47r2 



X2 + y2V ^X2 + F2 



Now, 



and 



Y{t - r) = 5r{t - r) = ri(t - r) - r{t - r) = ri(t) - r{t) = Y{t) = Y 



(C15) 



(C16) 



(C17) 
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X{t -t) = r{t - T)6e{t -t)= r{t - r) (Oiit - r) - e{t - t)) 



,,«,(()_«(,)_ I <l%(n)_(i%w,^ 



ri 



(C18) 



In the local approximation, we can expand the last term in equation (|C18|) in a Taylor series 
in ri — r. This yields 

X(t -t)= ridAt) - dit)) - r4- ff^^^}^] (n -r)r = X- J^MYt. (C19) 

dr \ r J 

where E(r) is the local shear of the flow ( p.l4| ). Substituting equations ( |C17] ) and ( |(J19D in 



equation 



we get 



C{r) = -ArP.Jr) I dXdY- 



X-i:{r)Yi 



47r2 ^ ""^^ ' J ^^^^^ X2 + Y^{X- ^{r)YTY + Y^' 
This integral is similar to (|C6|) , so we directly obtain equations (|113|) and (|118|) . 



(C20) 



APPENDIX D: CONSERVATION LAWS SATISFIED BY THE GENERALIZED 

KINETIC EQUATION 



In this section, we prove some general properties satisfied by equation (|128|) . Note first 
that it can be written 



dP 
'dt 



+ (V)VP = -V.J, 



where 



-A^ 



I dr I rfVV(l ^ 0)i|v(l ^ 0)PiVP + V(0 -^ l)PVPi| 



t-T 



(Dl) 



(D2) 



is the diffusion current. It is clear at first sight that equation ( |D 1| ) conserves the total circu- 
lation r = J{uj)d'^r provided that J.n = on the boundary. We now prove the conservation 
of other integral constraints depending on the domain shape. 

(i) In a circular and in an infinite domain, the angular momentum defined by 



{uj)r d r 



must be conserved. Taking the time derivative of equation 
remembering that L is conserved by the advective term, we get 



(D3) 
substituting for ( |nT| ) and 

L = 2iV7 / J.rd^r. (D4) 

Substituting explicitly for the diffusion current (|D2|) in equation (|D4| ) , we obtain 
L = -2N^-f I dr f dhdhi[v.Y{l -^ 0)]i|v(l -^ 0)PiVP + V(0 -^ l)PVPi| . (D5) 
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Permuting the dummy variables r and ri, we get 

L = -2N^-f I dr f dhdhi[ri.Y{0 -^ l)]t|v(l -^ 0)PiVP + V(0 -^ l)PVPi| . (D6) 
Adding these two quantities, we arrive at the final expression 

L 

x|v(1^0)PiVP + V(0^ l)PVPi| . (D7) 



-A^'7 / dr I d'rd'ri[r.V{l ^ 0) + ri.V(0 -^ l)]t 



Now, the term in brackets vanishes as shown by the following argument. Consider two point 
vortices in a circular (or infinite) domain. Their angular momentum is 



7(r^ + rl) 



and it is conserved. This implies: 



° - f - 2<4 + ^-^ ) - 27[r.V(l ^ 0) + r,V(0 -. 1)]. 



(D8) 



(D9) 



We can also prove this result by a direct calculation. In an unbounded domain V(0 — > 1) 
— V(l — » 0) and consequently 



r.V(l ^ 0) + ri.V(0 ^ 1) = V(l -^ 0).(r - n) 



0, 



(DIO) 



where we have used equation (|^ to get the last equality. In a circular domain, the velocity 
V(l -^ 0) is given by equation (|^) plus a term V;,(l -^ 0) which can be determined with 
the method of "images" . If R denotes the domain radius, we find 



7 



V.(l-0)^-.A 



Sri-r 



ri 



(Dll) 



Therefore, 



r.V,(1^0) + ri.V6(0^i: 



7 



zArij.r (zAr).ri 






I n 



R 



R r |2 



R' 



Noting that (z A ri).r = — (z A r).ri and that 



R n 

— ri r 

n R 



P' + 1^^' - 2r.ri 
R"^ 



we finally conclude that 



r.V(l -> 0) + ri.V(0 -^ 1) 



R r 
— r ri 

r R 



0. 



(D12) 



(D13) 



(D14) 



From this identity and from equation ( P7| ), it results that the kinetic equation (|128 
serves the angular momentum in a disk and in an infinite domain, i.e. L = 0. 



con- 



43 



(ii) In an infinite domain or in a channel, the hnear impulse 

P = /"r A {uj)7.dh (D15) 

must be conserved (in a channel extending in the x direction, only the component P^ of the 
linear impulse must be conserved). Taking the time derivative of equation ( pi5| ), substitut- 
ing for (pi|) and remembering that P is conserved by the advective term, we get 

P^ = N-f f Jdh. (D16) 



Substituting explicitly for the diffusion current (|D2|) in equation ( pl6| ), permutting the 
dummy variables r and ri, and taking the half-sum of the resulting expressions we finally 
obtain 



N'-f 



2o, rt 



2 

X 



f dr f dhdhi[V{l -^ 0) + V(0 -^ l)]t 
|v(1^0)PiVP + V(0^1)PVPi| . (D17) 



Now, we can use the same argument as before to show that the term in brackets vanishes. 
Let us consider two point vortices in a channel (or in an infinite domain). Their linear 
impulse is 

P± = 7(r + ri) (D18) 

and it is conserved. This implies 

dP± fdr dri" 







dt 



<i + iF)=^i^(^^'"+^(°-^)i- '^1^' 



We can also prove this result by a direct calculation. Equation ( P19| ) is obvious in an 



unbounded domain since V(0 ^ 1) = — V(l -^ 0). In a channel extending in the x 
direction, we need to show that Px = 0, i.e. Vy{l ^ 0) + Vy{0 ^ 1) = 0. Now, the velocity 
Vy{l -^ 0) is given by equation (|^) plus a term Vfe(l —* 0)y which can be determined with 
the method of "images". If a denotes the width of the channel, we find that 

^ zi2? r X — Xi X — Xi 1 

H(l - 0), = -— ^g^l ^^ _ ^^^2 + (y-y,- 2naf " (x - x{)^ + {y + yi- 2na)A' ^^^°^ 

Under this form, it is clear that 14(1 — » 0)^ is antisymmetric under the exchange of 1 and 
0, so that finally 

\4(1^0) + 14(0^1) = 0. (D21) 

From this identity and from equation (|D17|) , it results that the kinetic equation ( p.28| ) con- 
serves the linear impulse in a channel or in an infinite domain, i.e. P = 0. 

APPENDIX E: CALCULATION OF THE MEMORY FUNCTION 

In this Appendix, we calculate the memory function that occurs in equation ( p.28|) . 
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1. Axisymmetrical flow 

If we assume that P = P{r,t), then equation (|128D simphfies in 



dP Id 

where 

Jr = -N J^ drj d^riVr^t){l -^ 0)t[Vrit-r){l ^ 0)Pi— - V;,(t-,)(1 ^ 0)P^| (E2) 

and where K(j)(l — ^ 0) is the component of the vector V(l -^ 0) in the direction of r(t). If 
we denote by {r(t),9(t)) and {ri(t),9i(t)) the polar coordinates that specify the position of 
the point vortices and 1 at time t, we easily find that 



7 Ti sin{6 — 9i) 

2t[ r\ + r"^ — 2rri cos{6 — 6i) 






We shall assume that between t and t — r, the point vortices follow circular trajectories with 
angular velocity Vt{r,t). In that case, r{t — t) = r and 6{t — t) = 6 — fl{r,t)T. Then, we 
obtain 

^(*-"^^ ' 27rrl + r^-2rr,cos{9-ei-AnT) ^ ^ 

with 

An = n{r,t)-n{ri,t). (E5) 

We find similarly that Ki(t-r)(l ^ 0) = — K-(t_r)(l -^ 0). Our previous assumptions also 
imply that P{r{t — T),t — t) ~ P{r,t) between t and t — t. In words, this means that 
the correlation time is smaller than the time scale on which the average vorticity changes 
appreciably. We do not assume that it is much smaller, so this approximation is not over 
restrictive. In that case, the diffusion current becomes 



Jr = -N 



r-\-oo /'27r z'+oo 



r-\-oo rZTT r-\-oo 

/ dr dOi rridriKw(l-^0)K(t-r)(1^0) 

JO JO JO 



lp^9P _ ]_^m 

r dr ri dri 



(E6) 



where the time integral has been extended to +oo. We now need to evaluate the memory 
function 



M = J^ drj^ de^Vrit){l^O)Vrit-r){l^O)- (E7) 



Introducing the notations (p = 9i — 9 and 



rf + r 
we have explicitly: 
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2rr^ 



V47rr/ Jo Jo 1 - Acos0 1 - A cos(0 + Afir) ^ ^ 

This can also be written 

M=[^] / dr d(PV'U)V'U+^nT), (ElO) 

V47rr/ JO io 

where 

y(0) = ln(l-Acos0). (Ell) 

We now write the function V^(0) in the form of a Fourier series 

^(0) = E «ne'"^ with an = — \/(0)e-^"*#. (E12) 



The memory function becomes 

2 /■+00 /•27r 



M = — f^j / dr/ rf0 Y. rima„a^e*("+'")V"*^^". (E13) 



n,rra=— oo 



Carrying out the integrations on and r using the integral representation of the delta 
function 

5{x) = — / e-^^^dp, (E14) 

we are left with 



M = -— Y^ nmanamSn-mS{mAQ) = —S{AQ) ^ |ra|a^. (E15) 

n,m,=— oo ra=— oo 

It remains for us to evaluate the series that appears in the last expression of the memory 
function. Using the identities: 



rin(l-Acos0)cosM)rf0= --f-- W— -l) (n>0) (E16) 

Jo n\X V X'^ J 



ln(l - A cos <p)d<j) = Ti In (^ + — — — ] (E17) 



and the definition (|E8| ) of A, we find that a^ < oo and, for ri > 0: 

_ 1 /(rj + r^) - |rj -r^K" _ 1 /r<\'^ 



ra V 2rri J n Vr- ' 



(E18) 



where r> (resp. r<) is the biggest (resp. smallest) of r and ri. Therefore, the value of the 
series is 
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oo l^^^x2n 



n \r 



E \n\al = 2j:nal = 2j:^(:^ 

n=—oo n=l 

The memory function takes the form 



-2 In 



1- ^ 



(E19) 



M 



i-A^^)^^ 



1-1^ 



2i 



and the diffusion current in the axisymmetrical case can be written 



Jr = / ridri6(Q — Qi) In 

4r JO 

This leads to the kinetic equation ( |133[ ). 



1- 



r<\2] rl 5P 1 dPi 



-Pi 



P- 



r dr Ti dri 



(E20) 



(E21) 



2. Unidirectional flow 



If we assume that P = P{y,t), then equation (|128|) simphfies in 



dP 
'dt 



dy 



with 



J, = -N l\lr j ,l..,ly,V(l ^ 0),V'.(1 - 0)._,{/>,§^ - P^}^_. 



(E22) 



(E23) 



Assuming that between t and t — t the vortices follow linear trajectories with velocity 
(V) = {V){ii,t)ex, we have y{t — r) = y and x{t — r) = x — {V){y,t)T. Therefore, the 
function V^^(l — »• 0) at times t and t — r takes explicitly the form 



and 



vy{i ^ o)t 



V^il -> 0)t-r 



7 



Xi — X 



2rr {xi - xy + {yi - y^ 
7 Xi — X — AVt 



27r {xi-x- AVry + (t/i - y^ ' 
where we have introduced the notation 

AV={V)iy^,t)-{V)iy,t). 



(E24) 



(E25) 



(E26) 



We also assume that the correlation time is smaller than the time scale over which the 
vorticity changes appreciably. Then, P{y{t — T),t ~ t) ~ P{y,t) and the diffusion current 
becomes 



Jy = -N J^ drj dx,dy,Vy{l -^ 0)tVy{l -^ 0)i_.|Pi— - P^ 
We now need to calculate the memory function 



(E27) 
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n+OO 

M 



r+oo r+oo 

/ dr / dxiVy{l -^ 0)i1/^(l ^ 0)i_,. (E28) 

Using equations ([E24|) and (|E25|), we have explicitly 



7^ /•+00 .+00 X X - AVr 

where we have set X = xi — x, F = yi — y. Equation (|E29| ) can also be written 

with 



Vr(X,y) =lnVX2+y2 = ln^ (E31) 

We shall now write the function W{^ in the form of a Fourier integral 

W{^) = — ^ / W(k)e"^d\ with #(k) = / W{$,)e^^d% (E32) 

Then, the foregoing expression for the memory function becomes 

/ rfr / dX d''M^\i:W{\i)W{\i!)Kk',e-'^''+''^e'^^^'^\ (E33) 

JO J-oo J 



M -- ^ 



(27r)6 Jo 
Carrying out the integrations over r and X, we get 

„2 



M = --^ / d^M^\i!W{\i)W{\i!)k^k'J{k, + A;^.)e"'^''''+''^^^5(A;1.AK) (E34) 

327r'^ J 

and, consequently, 

/+00 
^ dk,dkydk'^W{k,, ky)W{-k,, k'y)\k4e-'^''y+''y^^ . (E35) 



M ' 



Now, the Fourier transform of W can be written explicitly 

f + CO 

W{k) = 27r / WiOMkO^d^, (E36) 

JO 

where use has been made of the well-known identity 

/ cos{zcose)de = 27rJo{z), (E37) 

JO 

where Jo is Bessel function of order zero. It is immediate to see that the Fourier transform of 
W as defined by (|E31|) does not exists. Indeed, the integral ( Pil36| ) diverges when ^ —>■ +cx3, 
i.e. at large separations. However, in physical situations the domain never extends to infinity 
so that, in practice, the integral remains finite. A convenient way to introduce a cut-off at 
large separations is to make the substitution 
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W(e)=lne - W{i) = -K,{i/L), (E38) 

where Kq is the modified Bessel function of order zero and L is a length scale of the order 
of the system size. For small separations Kq{z) ~ — Inz and for large separations Kq{z) ~ 
, /^e~^. This modification amounts to replacing the Poisson equation (Q) by an equation of 
the form 

- Ail) + -^ = uj (E39) 



Equation ( [E39D is precisely what is obtained in geophysics in the "quasigeostrophic approx- 



imation". The deformation of the fiuid surface introduces a shielding of the interaction 
between vortices on a length ~ L, called the Rossby radius. Obviously, the Rossby radius 
plays the same role as the Debye length in plasma physics. 
With this prescription we find that 

W^(k) = -^^ (E40) 



where ki = 1/ L. Substituting in equation ([E35|) , we get 



'' - 2^^'^^) IZ ''MC ^fT&*»)'- '"^^' 



The integral on ky can be carried out easily, leaving the result 
Now, setting t^ = AY'^{kl + k1), we finally obtain 



M = ^6{AV) r^ rffc,_^-e-2|^lv/^i+^. (E42) 

4 Jo kf + ki 



M=^5(A\/)Ei(^), (E43) 



where 



~oo g~* 



Ei{x) = / — dt (E44) 

Jx t 

is the exponential integral. For x — > 0, we have Ei{x) = — C — Inx where C = 0.57721... 
is Euler's constant. In conclusion, the diffusion current in the unidirectional case takes the 
form 

and it leads to the kinetic equation ( |135| ). 
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